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20.  cont. 

► reversibil  ity  and  chromatographic  effects  would  influence  the  design  and  subse- 
quent operation  of  prototype  adsorbers. 

Numerous  adsorption  equilibrium  studies  were  conducted  using  bi-solute  com- 
binations of  ^-phenyl phenol  (OPP),  dinitro-o-sec-butyphenol  (DNOSBP),  and 
2 ,4-dichlorophenol  (2,4-DCP).  Results  of  these  studies  were  employed  to  eval- 
uate the  utility  of  the  Longmuir  Competitive,  Longmuir  Semi-Competitive  Ideal 
Solution  Theory  and  Graphical  Models  of  competitive  adsorption  equilibria. 
Neither  of  the  models  was  found  to  adequately  describe  competitive  adsorption 
effects.  Nonetheless,  the  Graphical  approach  gave  the  best  estimate  of  com- 
petitive interactions.^. 
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Simulations  using  \he  dynamic  mathematical  model  showed  that: 

(1)  the  fluidizedJbed  mode  of  operation  offers  distinct  advantages  over 
packed-bed  operation  for  minimizing  the  elution  of  adsorbed  contami- 
nants from  adsorbers; 

(2)  industrial  production  schedules  should  be  arranged  such  that  all  con- 
taminants are  discharged  concurrently; 

(3)  production  schedules  should  be  arranged  such  that  wastewaters  con- 
taining the  highest  energy  adsorbing  contaminants  are  introduced  firs 
to  a column  of  fresh  adsorbent  while  the  wastewaters  containing  the 
lowest  energy  adsorbing  contaminants  are  introduced  last; 

(4)  concentration  equilization  should  be  employed  ahead  of  adsorbers  to 
minimize  displacement  and  subsequent  elution  of  adsorbed  contaminants 

Columnar  laboratory  studies  generally  verified  the  conclusions  drawn  from  the 
simulations. 
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SUMMARY 


As  water  quality  standards  have  become  more  stringent,  there 
has  been  an  increasing  interest  in  the  use  of  activated  carbon  for 
the  removal  of  trace  contaminants.  Accordingly,  the  objective  of 
this  study  was  to  develop  a predictive  model  for  adsorption  of  mul- 
tiple solutes  from  solution  onto  activated  carbon  using  a differen- 
tial contacting  system.  To  accomplish  this  goal,  two  secondary 
objectives  had  to  be  achieved.  These  were  (1)  to  evaluate  existing 
multi-solute  adsorption  equilibria  models  and  (2)  to  determine  how 
adsorption  reversibil ity  and  chromatographic  effects  would  influence 
and  subsequent  operation  of  prototype  adsorbers. 

Numerous  adsorption  equilibrium  studies  were  conducted  using 
bi-solute  combinations  of  o^-phenyl phenol  (OPP),  dinitro-o-sec-buty- 
phenol  (DNOSBP),  and  2 ,4-dichlorophenol  (2,4-DCP).  Results  of  these 
studies  were  employed  to  evaluate  the  utility  of  the  Longmuir  Com- 
petitive, Longmuir  Semi-Competitive  Ideal  Solution  Theory  and  Gra- 
phical Models  of  competitive  adsorption  equilibria.  Neither  of  the 
models  was  found  to  adequately  describe  competitive  adsorption  effects. 
Nonetheless,  the  Graphical  approach  gave  the  best  estimate  of  com- 
petitive interactions. 

Simulations  using  the  dynamic  mathematical  model  showed  that: 

(1)  the  fluidized-bed  mode  of  operation  offers  distinct  advan- 
tages over  packed-bed  operation  for  minimizing  the  elution 
of  adsorbed  contaminants  from  adsorbers; 

(2)  industrial  production  schedules  should  be  arranged  such 
that  all  contaminants  are  discharged  concurrently ; 

(3)  production  schedules  should  be  arranged  such  that  waste- 
waters  containing  the  highest  energy  adsorbing  contaminants 
are  introduced  first  to  a column  of  fresh  adsorbent  while  the 
wastewaters  containing  the  lowest  energy  adsorbing  contam- 
inants are  introduced  last; 

(4)  concentration  equilization  should  be  employed  ahead  of 
adsorbers  to  minimize  displacement  and  subsequent  elution 
of  adsorbed  contaminants. 

Columnar  laboratory  studies  generally  verified  the  conclusions  drawn 
from  the  simulations. 
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CHAPTER  I 


INTRODUCTION 

Use  of  adsorption  technology  for  the  decontamination  and  renovation 
of  wastewaters  has  during  the  past  decade  and  particularly  during  the 
past  five  years  been  implemented  for  full-scale  systems  with  increasing 
frequency.  This  is  because  adsorption  on  both  granular  and  powdered 
activated  carbon  has  clearly  been  shown  to  have  rather  unique  capabilities 
for  removal  of  trace  levels  of  pollutants,  such  as  pesticides  and  other 
surface  active  organic  materials,  from  wastewaters.  Throughout  the 
United  States  approximately  thirty  adsorption  installations  that  are  to 
be  employed  for  the  treatment  of  municipal  wastewaters  are  either  in  the 
desiqn  or  construction  stages  or  have,  indeed,  already  been  placed  on- 
line. Likewise,  numerous  industries  are  currently  using  adsorption 
contactors  either  for  processing  applications  or  for  removal  of  certain 
noxious  contaminants  from  its  aqueous  waste  stream. 

Continuous-flow  adsorption  contactors  have  heretofore  generally  been 
designed  for  operation  in  the  packed-bed  mode.  It  is  to  be  noted  that 
packed  beds  of  nranular  active  carbon  are  well  suited  for  the  treatment 
of  liquids  that  contain  little  or  no  suspended  solids  and,  with  a clear 
feed  stream,  can  be  expected  to  operate  effectively  for  extended  periods 
without  clogging  or  excessive  pressure  loss.  The  presence  of  suspended 
solids  in  aqueous  waste  streams,  however,  presents  serious  problems  for 
the  use  of  activated  carbon  in  packed  beds.  These  solids  lead  to  a 
progressive  clogging  of  the  beds,  much  as  in  a deep-bed  (sand)  filter, 
with  resulting  increases  in  head  loss.  Conversely,  a fluidized-bed 
adsorber,  in  which  the  wastewater  which  is  to  be  treated  is  passed  upward 
through  a bed  of  activated  carbon  at  a velocity  sufficient  to  expand 


the  bed,  is  not  plagued  with  problems  of  bed  clogging  and  increasing 
pressure  loss.  Whereas,  the  former  contacting  scheme  serves  in  a dual 
capacity  as  a filter  as  well  as  an  adsorber,  the  latter  functions  only 
in  the  singular  role  as  an  adsorption  process. 

Although  the  performance  of  packed-  and  fluidized-bed  adsorption 
contactors  has  been  observed  to  be  virtually  identical  when  employed 
for  the  treatment  of  municipal  waste  waters,  substantial  differences  in 
column  performance  can  be  anticipated  for  the  two  contacting  modes  when 
industrial  waste  waters  having  a well  defined  composition  of  high  and 
low  energy  adsorbing  solutes  are  treated  by  the  adsorption  process. 

It  is  clear,  therefore,  that  the  choice  between  the  two  contacting 
systems  is  heavily  deDendent  upon  the  specific  solution  and  particulate 
characteristics  of  the  wastewater  that  is  to  be  treated.  The  choice 
cannot  be  defined  by  a cursory  analysis,  but  can  only  be  predicted  on  a 
detailed  investigation  of  the  characteristics  of  the  waste  stream  that 
is  to  be  treated  and  by  subsequent  simulation  of  the  anticipated 
column  performance. 

Both  design  and  operation  of  adsorption  contactors  are  related  to 
the  characteri sti c breakthrough  curves  in  which  the  readily  measured 
effluent  concentration  is  plotted  as  a function  of  elapsed  time. 
Depending  on  the  ultimate  objective  of  the  process,  either  the  break- 
throunh  Doint  or  the  saturation  point  represents  the  design  and/or 
operational  target.  Prediction  of  the  breakthrough  phenomenon  has 
occupied  the  attention  of  numerous  investiaators  during  recent  years. 
Nonetheless,  throughout  this  period  specification  of  design  and 
operational  criteria  for  full-scale  adsorbers  has  been  predicated 
principally  on  empirical  observations.  This  is  primarily  due  to  the 
fact  that  municipal  wastewaters  have  an  extremely  complex  composition 
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which,  in  turn,  substantially  complicates  modeling  and  simulation  pro- 
cedures. For  this  situation  any  practicable  model  must  of  necessity 
be  of  a quasi-empi rical  nature  to  be  generally  applicable. 
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In  contrast,  the  composition  of  industrial  wastewaters  frequently 
can  be  explicitly  defined.  For  example,  the  wastewaters  resulting 
from  the  manufacture  or  handling  of  pesticides  would  likely  consist  of 
several  major  constituents.  While  the  concentration  of  each  such  con- 
taminant would  undoubtedly  be  time  variant,  the  composition  could 
nevertheless  be  clearly  determined.  In  this  case  a generalized  descriptive 
model  can  be  employed  for  accurately  simulating  the  dynamic  responses  of 
columnar  adsorbers  for  a variety  of  operational  and  systemic  conditions 
including  adsorber  geometry  and  configuration,  flow  rate,  composition, 
concentration,  temperature,  contacting  modes,  competitive  effects, 
reversibi 1 i ty , etc.  Use  of  such  a model  which  is  wholly  fundamental  in 
scone  provides  one  with  a rational  basis  for  evaluating  design  and 
operational  al ternati ves , for  projecting  optimal  ranges  of  design  and 
operational  criteria,  and  for  simulating  process  control  strategies. 
Moreover,  such  a mathematical  development  is  extremely  beneficial  for 
SDecification  of  the  cost  effective  design  of  a pilot-scale  system  if 
such  were  indicated.  Under  certain  circumstances,  however,  use  of  a 
fully  descriptive  model  should,  indeed,  be  entirely  satisfactory  for 
specifying  design  and  operational  parameters  without  the  need  to  resort 


to  time-consuming  and  costly  pilot-scale  studies. 
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CHAPTER  II 

OBJECTIVES  OF  RESEARCH 


Numerous  research  projects  have  demonstrated  the  feasibility  of 
utilizing  the  unit  process  of  adsorption  on  activated  carbon  for  the 
removal  of  pesticides  from  waters.  Others  have  delineated  and 
characterized  process  parameters  significant  to  the  performance  of 
packed-bed  and  fluid-bed  activated  carbon  adsorbers.  While  these 
investigations  have  demonstrated  the  unique  utility  of  the  process 
and  have  yielded  useful  information  regarding  general  design  and 
operating  criteria  for  pesticide  decontamination  applications,  there 
has  remained,  however,  a definite  need  for  the  systematization  of  this 
information,  and  for  the  development  of  a quantitative  rational  design 
procedure  for  adsorbers  for  field  applications. 

Because  prediction  or  simulation  of  tne  concentration-time  profiles 
that  characterize  adsorber  performance  is  paramount  to  the  development 
of  rational  design  and  operational  criteria,  various  theoretical  math- 
ematical models  have  previously  been  developed  for  predicting  breakthrough 

profiles  for  continuous-flow  adsorbers  in  which  sinqle-solute  solutions 
were  contacted  with  the  solid-phase.  Subsequently,  a study  initiated 
and  directed  by  the  principal  investigator  resulted  in  the  development 
of  a mathematical  model  which  can  be  employed  to  simulate  the  performance 
of  either  packed-  or  fl uidi zed-bed  adsorption  contactors  when  contacted 
with  multiple-solute  mixtures  such  as  are  commonly  encountered  in  indus- 
trial wastewater  discharges.  This  model, in  part,  accounted  for  the 
competitive  effects  between  solutes  with  regard  to  adsorption  equilibria. 
Further,  according  to  fundamental  thermodynamic  principles  the  model  was 
formulated  on  the  assumption  that  adsorption  on  activated  carbon  was 
fully  reversible.  This,  of  course,  had  been  postulated  by  numerous 
previous  investigators. 
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Experimental  investigations  which  were  conducted  durinq  the  latter 
portion  of  the  contract  period,  however,  showed  that  the  assumption  of 
full  reversibility  was  not  entirely  valid  for  description  of  adsorber 
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column  performance  in  the  dynamic  state.  The  results  showed  neither 
total  reversibility  nor  total  irreversibility  in  the  time-scale  that  is 
significant  to  adsorber  operation;  although  it  is  anticipated  that  if 
an  adsorption  column  were  allowed  to  equilibrate  over  very  extended 
periods  of  time  one  would,  indeed,  observe  a fully  reversible  condition. 
While  the  latter  condition  is  significant  with  regard  to  adsorption 
thermodynamics,  it  has  little  practical  significance  regarding  description 
of  the  performance  of  adsorption  contactors  which  would  be  employed  for 
the  removal  of  residues,  such  as  pesticides,  from  wastewaters.  In  this 
case  any  fully  descriptive  model  should  account  for  the  relative  extent 
of  reversibility/irreversibility  in  the  time-scale  that  has  engineering 
significance.  The  essence  of  this  research,  therefore,  is  to  sufficiently 
modify  the  existing  mathematical  developments  to  account  for  all  time- 
dependent  reversibi 1 i ty/i rreversibi 1 i ty  considerations  such  that  mathema- 
tical simulations  of  adsorber  column  performance  accurately  predict  these 
time  dependent  phenomena.  Accordingly,  the  specific  aims  and  objectives 
of  this  investigation  are  enumerated  as  follows: 

1)  To  modify  an  existing  theoretical  mathematical  model,  which  was 
formulated  under  the  direction  of  the  principal  investigator  to 
accommodate  the  effects  of  time-dependent  relative  reversibility/ 
i rreversibi 1 i ty ; 

2)  To  conduct  batch-system  equilibrium  studies  for  single-solute 
systems  using  pesticide-type  compounds  that  are  of  direct 
interest  to  the  Department  of  the  Army.  Single  solute  adsorption 
parameters  so  obtained  can  then  be  directly  employed  for  pre- 
diction of  multiple-solute  adsorption  equilibria  which  includes 
competitive  effects; 
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3)  To  experimentally  determine  the  time-dependent  reversible/ 
irreversible  interactions  between  solutes  in  mul ti -component 
systems.  These  experiments  were  to  be  conducted  in  columnar 
continuous-flow  contactors  so  as  to  account  for  the  time- 
dependent  characteristics  of  reversibility; 

4)  To  simulate  the  multiple-solute  adsorption  process  through 

use  of  the  model  formulated  under  Part  1.  Experimental  verification 
is  to  be  obtained  from  a series  of  columnar  adsorption  studies; 

5)  Upon  obtaining  satisfactory  agreement  between  the  experimental 
results  and  model  simulations,  to  conduct  an  extensive  series 

of  simulations  to  fully  evaluate  all  possible  dynamic  conditions 
that  might  likely  be  encountered  in  field  applications.  This 
serves  a two-fold  purpose:  (a)  to  define  which  design  and 

operational  alternatives  would  provide  for  optimal  performance 
of  activated  carbon  adsorbers  which  are  subjected  to  specific 
and  critical  dynamic  input  conditions;  and  (b)  to  illustrate 
how  one  would  proceed  to  employ  the  mathematical  developments 
to  simulate  design,  operational,  and  process  control  alternatives. 

Upon  completion  of  the  objectives  detailed  above.  Department  of  Army 
scientists  and  engineers  will  have  at  their  disposal  a relatively  straight- 
forward method  for  design  of  adsorption  contactors,  for  specifying  opera- 
tional criteria  for  such  adsorbers,  and  for  defining  process  control 
strategies.  This  will  under  most  circumstances  obviate  the  need  for  lengthy 
and  costly  pilot-scale  studies  which  heretofore  have  been  reguired  before 
design  could  be  initiated.  Alternatively,  if  because  of  special  conditions, 
pilot  studies  are  indeed  indicated,  then  the  model  developments  can  be  employed 
to  plan  the  experimental  pilot  studies  so  as  to  maximize  the  quantity  of 
useful  data  obtained  per  unit  of  cash  expended.  Both  aspects,  obviously, 


will  serve  to  minimize  the  expenditures  required  of  the  Department  of 
Army  to  arrive  at  a cost  effective  and  fail-safe  design  of  adsorption 
contactors  for  specific  appl ications. 
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‘CHAPTER  III 
EXPERIMENTAL  APPROACH 


The  proposed  research  project  was  divided  into  five  phases.  These 
are  detailed  as  follows: 

PHASE  I:  MODEL  MODIFICATION  AND  REFINEMENT 

Under  a previous  contract  (43)  a mathematical  model  was  developed  that 
is  general  and  applicable  for  the  prediction  of  the  performance  of  adsorbers 
when  treating  wastewaters  that  contained  several  major  contaminants.  The 
model  was  based  on  general  mass  balance  relationships  and  the  appropriate 
rate  and  equilibrium  parameters.  With  regard  to  adsorption  kinetics  the 
model  was  formulated  assuming  that  film  diffusion  was  the  predominant  rate 
limiting  mechanism.  This  assumption  is  well  supported  in  the  literature 
(see  Chapter  5).  If,  in  contrast,  pore-diffusion  or  pore-surface  diffusion 
were  determined  to  be  rate  limiting  for  a specific  application,  the  model 
was  sufficiently  general  to  accommodate  these  alternate  limiting  mechanisms. 
Further,  competitive  adsorption  parameters  that  were  included  in  the  model 
were  based  on  fundamental  thermodynamic  principles  which  were  found  to  ac- 
curately predict  the  competitive  effects  of  adsorption  equilibria. 

One  other  assumption  that  was  employed  during  model  formulation  was 
that  adsorption  on  granular  activated  carbon  was  fully  reversible.  Again, 
this  assumption  was  extensively  supported  in  the  literature,  particularly 
in  that  literature  which  pertains  to  chromatography.  Nevertheless,  several 
mathematical  formulations  which  had  been  reported  in  the  literature  were 
based  on  the  assumption  of  complete  irreversibil ity.  This  latter  assumption, 
however,  was  primarily  for  the  purpose  of  mathematical  expediency  such 
that  an  exact  solution  of  the  set  of  equations  could  be  effected.  As  was 
indicated  previously,  experimental  investigations  that  were  conducted  under 
the  previous  contract  clearly  showed  that,  in  the  dynamic  state  and  in  the 
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time  domain  of  significance  in  the  operation  of  adsorption  columns, 
adsorption  of  a variety  of  pesticides  was  neither  fully  reversible 
nor  completely  irreversible.  Accordingly,  the  principal  focus  of 
this  phase  of  the  investigation  was  to  incorporate  the  reversibility/ 
i rreversibi 1 i ty  considerations  into  the  model  which  was  previously 
developed. 

PHASE  II:  ACQUISITION  OF  COMPETITIVE  ADSORPTION  PARAMETERS 

Because  a theoretical  interpretation  of  multicomponent,  competitive 
adsorption  equilibria  was  obtained  during  the  previous  contract,  the 
need  to  perform  additional  multiple-solute  equilibrium  experiments  was 
obviated.  Consequently,  only  single-solute  experiments  were  conducted. 
These  provided  data  that  is  suitable  for  prediction  of  the  competitive 
effects  of  mul ti-component  adsorption  equilibria.  These  were  conducted 
for  three  organic  herbicides;  dinitro-o-sec-butylphenol  (DNOSBP), 
o-phenylphenol  (OPP),  and  2 ,4-dichlorophenol  (2,4-DCP). 

All  sinqle-solute  adsorption  data  were  obtained  using  four  con- 
tinuous stirred-tank  reactors  (CSTR)  operated  in  parallel.  Inlet 
concentrations  to  the  reactors  were  controlled  at  predetermined  levels. 

In  this  manner,  the  inlet  concentration  then  becomes  the  equilibrium 
concentration  for  the  particular  solid-solvent-solute  system.  Equilibrium 
parameters  obtained  from  this  apparatus  were  incorporated  in  the  basic 
mathematical  model  for  simulation  of  effluent-time  profiles  for  columnar 
flow-through  systems. 
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PHASE  III:  EXPERIMENTAL  DETERMINATION  OF  THE  PE  VERS IB I LI TY  OF 

COMPETITIVE  ADSORPTION 

This  phase  of  the  investigation  was  conducted  to  determine  the  time- 
dependent  reversible/irreversible  interactions  between  solutes  in  multi- 
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component  systems.  The  solutes  selected  in  Phase  II  of  the  study  were 
also  employed  in  this  facet  of  the  investigation.  Columnar,  continuous- 
flow  contactors  which  were  constructed  under  a previous  contract  were  used  to 
elucidate  the  parameters  that  are  required  to  describe  the  reversibility 
phenomena. 

With  one  exception,  the  experimental  procedure  involved  first 
saturating  the  carbon  with  a single  solute  in  the  columnar,  continuous- 
flow  contactors  emDloyed.  When  the  contactor  attained  a steady-state 
operational  level  (i.e.,  the  effluent  concentration  of  the  solute  was 
equal  to  the  influent  concentration  of  that  solute),  a second  solute 
at  the  desired  concentration  level  was  introduced  to  the  column  along 
with  the  first  solute  (concentration  was  not  changed  during  the  step 
input  forcing).  The  effect  of  the  second  solute  on  the  relative  column 
effluent  concentrations  was  then  observed.  The  run  was  continued 
until  steady-state  was  again  attained. 

In  one  experimental  study,  however,  both  solutes  (DNOSBP  and  OPP) 
were  applied  as  single  solutes  in  consecutive  order.  That  is,  DNOSBP 
was  applied  to  the  column  initially  until  steady-state  operation  was 
obtained.  OPP  was  subsequently  fed  to  the  column  as  a pure  solute. 

The  study  was  continued  until  steady-state  operation  was  again  attained. 

An  executive  synopsis  of  the  experimental  studies  conducted  follows: 

Fluidized-Bed  Contactor  Studies 

(a)  Application  of  both  DNOSBP  and  2,4-DCP  at  a molar  ratio  of 
1.0:0.913  after  presaturation  with  DNOSBP  (79.56  pm). 

(b)  Application  of  both  DNOSBP  and  2,^-DCP  at  a molar  ratio  of 
1.0:1.53  after  presaturation  with  DNOSBP  (80.0  pm). 

(c)  Application  of  both  2,4-DCP  and  DNOSBP  at  a molar  ratio  of 
1.0:1.007  after  presaturation  with  2,4-DCP  (80.0  urn). 
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(d)  Application  of  both  2,4-DCP  and  DNOSBP  at  a molar  ratio  of 
1.0:2.086  after  presaturation  with  2,4-DCP  (80.0  pm). 
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(e) 

(f) 

(g) 

(h) 

(i) 

(j) 

(k) 


Application  of  OPP  at  a concentration  of  80.0 
presaturation  with  DNOSBP  (80.0  pm). 
Application  of  both  OPP  and  DNOSBP  at  a molar 
after  presaturation  with  DNOSBP  (80.0  pm). 
Application  of  both  OPP  and  DNOSBP  at  a molar 
after  presaturation  with  DNOSBP  (80.0  pm). 
Application  of  both  OPP  and  DNOSBP  at  a molar 
after  presaturation  with  DNOSBP  (80.0  pm). 
Application  of  both  DNOSBP  and  OPP  at  a molar 
after  presaturation  with  OPP  (80.0  pm). 
Application  of  both  DNOSBP  and  OPP  at  a molar 
after  presaturation  with  OPP  (80.0  pm). 
Application  of  both  DNOSBP  and  OPP  at  a molar 
after  presaturation  with  OPP  (80. 0 pm). 


pm  after 


ratio  of  1.0: 2.0 


ratio  of  1. 0:1.0 


rati o of  1.8:1 


ratio  of  1. 0:2.0 


ratio  of  1.0: 1.0 


rati o of  1. 8: 1. 0 


Packed-Bed  Contactor  Studies 

(a)  Application  of  both  DNOSBP  and  2,4-DCP  at  a molar  ratio  of  1.0: 
0.42  after  presaturation  with  DNOSBP  (79.88  pm) 

(b)  Application  of  both  DNOSBP  and  2,4-DCP  at  a molar  ratio  of 
1.0:1.56  after  presaturation  with  DNOSBP  (78.52  pm). 

(c)  Application  of  both  2,4-DCP  and  DNOSBP  at  a molar  ratio  of 
1. 0:1.0  after  presaturation  with  2,4-DCP  (80.0  pm). 

(d)  Application  of  both  2,4-DCP  and  DNOSBP  at  a molar  ratio  of 
1. 0:2.0  after  presaturation  with  2,4-DCP  (80.0  pm). 
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(a)  Application  of  both  DNOSBP  and  2,4-DCP  at  a molar  ratio  of 
1.0:0.96  after  presaturation  with  DNOSBP  (80.0  pm). 

(b)  Application  of  both  DNOSBP  and  2,4-DCP  at  a molar  ratio  of 
1.0:2.02  after  presaturati on  with  DNOSBP  (30.0  pm). 

(c)  Application  of  both  2,4-DCP  and  DNOSBP  at  a molar  ratio  of 
1.0:0.96  after  presaturation  with  2,4-DCP  (80.0  pm). 

(d)  Application  of  both  2,4-DCP  and  DNOSBP  at  a molar  ratio  of 
1.0:2.03  after  presaturation  with  2,4-DCP. 

PHASE  IV:  SIMULATIONS  ASP  '.'ERIE I CATIONS  OF  THE  MATHEMATICAL  MODEL 

Upon  completion  of  Phases  I,  II,  and  III  the  mathematical  model  as 
modified  in  Phase  I was  utilized  to  simulate  columnar  continuous-flow 
experiments.  Once  satisfactory  agreement  between  experimental  results 
and  model  simulations  was  obtained,  a series  of  simulations  was  conducted. 
These  serve  to  fully  evaluate  possible  dynamic  conditions  that  might 
likely  be  encountered  in  field  applications.  Certainly,  this  provides 
guidance  for  the  ultimate  users  of  the  mathematical  developments  and 
pointed  out  design  and  operational  targets. 
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CHAPTER  IV: 

EXPERIMENTAL  APPARATUS  AND  TECHNIQUES 


rHY^::AL  apparatus  ahd  operating  procedure 

All  experimental  runs  were  conducted  either  in  fluidized-bed  contactors, 
packed-bed  contactors  and  continuous  stirred  tank  reactors  (CSTR).  All  three 
modes  of  contact  were  used  in  order  to  make  possible  model  comparisons  for 
the  various  systems. 

FhOKED  and  fluidized  bed  systems 

Physically,  this  experimental  apparatus  consisted  of  four  pyrex  glass 
columns,  two-inches  in  diameter  and  five  feet  in  length.  Each  contained  150 
grams  of  the  granular  activated  carbon.  Each  column  is  fitted  at  its  base 
with  a distribution  chamber  filled  with  glass  beads  which  serves  to  dissipate 
the  velocity  head  at  the  influent  section,  thus  ensuring  steady-state,  one- 
dimensional flow  along  the  length  of  the  column.  A process  flow  diagram 
of  one  column  is  shown  in  Figure  IV  - 1 . 

The  make-up  water  system  consisted  of  a 75  gallon  make-up  water  storage 
tank,  a Chemtrix  Incorporated , Model  45  pH  controller,  and  three  variable 
speed  Masterflex  Model  7545-13  chemical  feed  pumps,  two  of  which  were  used 
in  conjunction  with  the  pH  controller.  These  pumps, depending  upon  the  pH 
signal,  fed  0.02M  hydrochloric  acid  and  0.02M  sodium  hydroxide  into  the  influent 
line  of  the  make-up  water  tank  at  rate  of  10  milliliters  per  minute  to  main- 
tain the  set-point  pH.  The  other  chemical  feed  pump  which  had  the  same 
capacity  as  the  other  pumps  was  used  to  feed  a phosphate  buffer  to  maintain 
a concentration  of  0.002  M.  The  make-up  water  was  delivered  to  each  of  the 
mixing  tanks  by  gravity  through  a one-half  inch  diameter  PVC  pipe  manifold 
system. 
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Flow  of  make-up  water  into  the  mixing  tank  was  controlled  by  a 3/8- 
inch  by  1/8-inch  solenoid  valve.  Operation  of  this  valve  was  controlled 
by  a level  control  system.  This  system  consisted  of  three  electrodes  in 
each  tank,  a low  and  high  water  elevation  electrode,  and  a ground. 

The  sorbate  feed  system  consisted  of  four  variable  speed  Masterflex 
pumps  (Model  7545-13)  which  were  capable  of  delivering  from  1.8  to  36 
milliliters  per  minute  or  0.06  milliliters  per  pump  revolution.  Tachometers 
were  attached  to  each  pump  to  monitor  speed  which,  in  turn,  was  related  to 
volumetric  flow  rate.  Normal  operational  speeds  for  these  pumps  were  120  RPM. 
This  is  equivalent  to  7.2  milliliters  per  minute.  Stock  solutions  were  supplied 
to  the  chemical  feed  pump  at  concentrations  of  either  2500,  5000  or  10,000 
micromoles/1 i ter. 

Solution  was  fed  to  the  columns  by  two  variable  speed  Masterflex  pumps 
Model  7545-17,  each  of  which  was  equipped  with  dual  pump-heads  capable  of 
delivering  84  to  1680  milliliters  per  minute.  This  capacity  was  such  that 
full  fluidization  of  an  activated  carbon  bed  of  400  grams  could  be  obtained. 

For  all  experimental  runs  the  flow  rate  through  the  contactor  was  maintained 
constant  at  5 gal lons/ft^/min  (0.413  liters/minute). 

The  sample  collection  system  consisted  of  a 30  gallon  receiving  reservoir, 
a solenoid  valve,  and  a cam  timer. 

Composite  samples  were  collected  on  the  basis  of  65  milliliters  per 
minute  for  a twenty-four  hour  period.  The  collection  vessels  were  then 
drained  at  the  end  of  sampling  periods  and  sampling  again  initiated.  Discrete 
samples  were  collected  at  twelve  hour  intervals.  Since  fluctuations  in  effluent 
were  gradual,  there  was  no  need  for  more  frequent  sampling. 
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The  packed  and  fluidized-bed  modes  of  operation  were  identical  with 
the  exception  of  the  direction  of  flow  through  the  column.  For  the  packed- 
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bed  mode  the  two  three-way  valves  were  positioned  such  that  the  flow  entered 
at  the  top  and  exited  through  the  bottom  of  the  column.  The  valve  positions 
were  reversed  for  the  fluidized-bed  mode  of  operation. 

CONTINUOUS  STIRRED  TANK  REACTOR  SYSTEM 

The  apparatus  used  consisted  of  four  nine-liter  baffled-reactors.  Figure 
IV-2  shows  the  process  flow  sheet.  The  feed  solution  was  fed  by  two  variable 
speed  Masterflex  pumps  Model  7545-15  with  dual  pumping  heads  each  of  which 
had  a rated  capacity  of  30  to  860  milliliters  per  minute.  At  the  onset 
of  the  runs  with  the  CSTR's,  a flow  rate  of  0.7  liters  per  hour  was  selected. 
Later  in  the  study  the  rate  was  increased  to  approximately  1.3  liters  per  hour 
due  to  the  extremely  long  periods  of  time  required  to  reach  saturation  at 
the  lower  flows.  Some  difficulty  was  experienced  in  maintaining  a steady  flow 
rate  at  these  low  levels.  For  that  reason  the  total  volume  of  effluent 
collected  from  the  reactor  was  used  to  compute  the  flow  rate  and  complete 
the  solid-phase  concentration  calcul ations  . 

The  rotameters  were  Gilmont  Model  3204-2  which  were  calibrated  for  flows 
of  10  to  860  milliliters  per  minute.  The  calibration  curve  for  each  of  these 
serial  numbered  rotameters  was  checked  with  a NBS  calibrated  Brooke  rotameter. 
All  the  calibration  curves  were  found  to  be  accurate. 

All  fittings  and  appurtenances  used  on  the  apparatus  were  brass,  produced 
by  Imperial  Eastman  Corporation.  All  tubing  was  3 8 inch  diameter 
polyethylene  tubing  (poly-flo),  manufactured  by  Imperial  Eastman  Corporation 
which  met  ASTM  0-1248  Type  1 Class  A,  Grade  4 specifications.  The  mixing 
system  for  each  reactor  consisted  of  a one-sixth  horsepower  motor,  stain- 
less steel  drive  shaft,  and  two  three-inch  diameter  impellers.  The  impeller 
consisted  of  six  blades  having  a pitch  of  45  degrees  with  the  vertical. 

The  effluent  system  consisted  of  a cylindrical  screen  3/8  inches 


in  diameter  and  approximately  4 inches  in  length.  It  was  constructed  of  a 
base  of  100  mesh  brass  alloy  wire  cloth  to  provide  structural  rigidity. 

A covering  of  325-mesh  brass  alloy  wire  cloth  was  placed  over  the  100  mesh 
material  to  ensure  retention  of  the  granular  activated  carbon. 

To  operate  this  system  on  a continuous  basis,  sixty  gallons  of  feed 
solution  had  to  be  prepared  prior  to  experimental  work.  The  feed  solution 
was  prepared  either  from  2500,  5000,  or  10,000  micromolar  stock  solution. 

The  pH  of  the  feed  solution  was  then  adjusted  to  the  pH  set-point  with  IN 
HC1 . The  system  was  then  buffered  with  a 0.002M  phosphate  buffer.  Samples 
of  this  feed  solution  were  analyzed  for  proper  concentration  and  pH.  The 
feed  solution  was  then  transferred  to  the  five  gallon  influent  tanks.  The 
solution  was  pumped  into  the  reactor  which  contained  10  grams  of  20/40  mesh 
activated  carbon.  The  mixers  were  operated  at  1000  RPM  which  ensured  that 
all  particles  remained  suspended  in  the  solution. 

Both  grab  samples  and  composite  samples  were  collected.  The  total 
effluent,  which  was  the  composite  sample,  was  collected,  measured  volumetri- 
cally,  and  the  concentration  determined.  The  influent  and  effluent  volu- 
metric measures  did  not  vary  more  than  one  percent.  This  indicates  the  close 
degree  of  control  obtained.  As  equilibrium  was  approached,  frequent  grab 
samples  were  taken.  When  the  concentration  in  influent  and  effluent  grab 
samples  were  equal , the  feed  solution  was  turned  off  and  the  system  was  allowed 
to  stir  for  twenty-four  hours.  If  the  concentration  in  the  reactor  did  not 
change  more  than  the  accuracy  of  the  analytical  procedure,  the  system  was 
considered  to  be  at  equilibrium.  The  composite  samples  were  then  used  to 
determine  the  quantity  of  sorbate  adsorbed  at  equilibrium. 

ADOOHBENT 

Filtrasorb  200,  manufactured  by  Calgon  Corporation,  was  the  activated 
carbon  used  in  all  studies.  Specifications  given  by  the  manufacturer  have 


been  detailed  in  Appendix  A. 

Prior  to  using  the  carbon  it  was  necessary  to  sieve,  wash,  and  dry 
the  carbon.  Since  the  manufacturer  provided  a wide  range  of  particle  sizes, 
it  was  necessary  to  select  a specific  particle  size  range  for  this  study. 

It  was  decided  that  a sieve  cut  of  that  passing  the  20  mesh  sieve  and  re- 
tained on  the  40  mesh  sieve  would  be  appropriate  for  the  studies. 

After  the  carbon  was  sieved  and  the  cuts  selected,  the  carbon  was 
washed  three  times  in  distilled  water  to  remove  the  fines.  It  was  then  dried 
at  105°C  for  twenty-four  hours.  After  allowing  the  carbon  to  cool,  it  was 
stored  in  glass  containers  at  37°C. 

ADSORBATES 

The  solutes  2, 4-di ni tro-o-sec-butyl phenol  (DNOSBP),  ortho-phenyl - 
phenol  (OPP),  and  2,4-dichlorophenol  (2,4-DCP)  were  used  for  this  study. 

DNOSBP  was  selected  for  use  because  of  the  extensive  work  done  with  it  and 
its  well  documented  properties.  OPP  and  2,4-DCP  were  selected  for  several 
reasons  including  molecular  size  as  compared  to  DNOSBP  and  the  value  of  the 
acid-base  equilibrium  constant  as  compared  to  DNOSBP,  This  selection  provided 
compounds  with  markedly  different  structures  and  charge  characteristics. 

This,  of  course,  would  affect  the  competition  of  the  solutes  for  adsorption 
sites  on  the  activated  carbon. 

The  DNOSBP  was  obtained  in  commercial  grade  and,  therefore,  required 
recrystallization  for  purification.  This  was  accomplished  by  heating 
the  compound  to  the  melting  point  in  a hot-water  bath  and  then  adding  an 
equal  volume  of  95  percent  ethanol.  The  mixture  was  heated  to  boiling. 
Distilled  water  was  then  added  until  precipitation  of  the  material  was  ob- 
served. The  solution  was  then  placed  in  a freezer  and  was  stirred  frequently. 
As  the  sorbate  precipitated  from  the  ethanol -water  mixture,  the  solution 
was  decanted.  Upon  total  crystallization  the  precipitate  was  removed 
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by  filtration  and  dried  under  vacuum.  The  melting  point  was  determined 
for  the  DNOSBP  which  had  twice  been  recrystallized.  The  melting  point 
curve  for  the  recrystallized  DNOSBP  had  a tolerance  of  + 1°C.  The  OPP 
and  2,4-DCP  were  acquired  in  reagent  grade  and  required  no  purification. 

All  three  solutes  are  relatively  insoluble  and,  therefore,  required  high 
temperatures  and  pH's  to  render  them  soluble.  Sodium  hydroxide  was  employed 
to  elevate  the  pH. 

ANALYTICAL  TECHNIQUES 

All  solutions  and  samples  were  analyzed  spectrophotometri cal ly  using 
a Beckman  Acta,  Model  III,  spectrophotometer.  To  determine  the  concentra- 
tions of  the  comp 'unds  considered  it  was  necessary  to  establish  the  maximum 
adsorption  wavelengths  of  the  solutes.  DNOSBP  in  the  anion  form  has  a 
maximum  absorbance  at  a wavelength  of  376  nanometers.  Similarly,  the 
maximum  absorbance  wavelength  for  neutral  2,4-DCP  was  determined  to  be 
285  nanometers  and  that  for  neutral  OPP  was  found  to  be  282  nanometers. 

The  compounds  were  studied  as  bisolute  pairs  of  DNOSBP  and  OPP 
and  DNOSBP  & 2,4-DCP  at  system  pH's  of  7.2  and  6.1,  respectively,  One 
of  the  compounds,  DNOSBP,  exists  in  the  anion  form  and  the  other,  2,4- 
DCP  or  OPP,  exists  in  the  neutral  form  at  the  pH's  studied. 

Analysis  of  the  DN0SBP/2,4-DCP  system  was  accomplished  through  the 
application  of  Beer's  Law  as  it  applies  to  two  solutes.  Accordingly,  it 
was  necessary  to  determine  the  extinction  coefficients  for  each  compound 
at  the  other  competing  solute's  maximum  absorbance  wavelength.  Table  I V -1 
lists  the  applicable  spectrophotometri c characteri sts  of  DNOSBP,  OPP,  and 
2,4-DCP. 

Analysis  of  the  DN0SBP/0PP  system  was  achieved  as  follows:  (1)  add  a 


volume  of  phosphate  buffer  to  one  portion  of  the  sample  to  poise  the  pH  at 
7.2;  (2)  add  a volume  of  sodium  hydroxide  to  another  portion  of  the  sample 
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to  raise  the  pH  to  11.2;  (3)  obtain  the  absorbance  reading  for  the  sample 
at  pH  = 7.2  at  a wavelength  of  376  nanometers  using  a distilltd  water  blank, 
and  (4)  obtain  the  absorbance  reading  for  the  sample  at  pH  = 11.2  at  a wave- 
length of  310  nanometers  using  the  sample  at  pH  = 7.2  as  the  blank.  The 
DNOSBP  concentration  was  determined  from  the  absorbance  reading  obtained 
at  pH  = 7.2  while  the  OPP  concentration  was  calculated  from  the  adsorbance 
obtained  at  pH  = 11.2. 
icedure 

Samples  of  approximately  100  milliliters  were  taken  from  the  influent, 
effluent,  and  sample  collection  tank  of  the  system.  Influent  and  effluent 
grab  samples  were  taken  during  a run  with  the  frequency  of  sampling  being 
greater  at  the  outset  of  the  run.  Samples  of  the  composite  were  taken  at 
the  beginning  of  each  day.  These  were  used  to  keep  a "running"  check  on  the 
olid-phase  concentration  profiles  by  having  available  a 24-hour  sample 
representative  of  the  average  effluent  for  that  day.  The  analysis  for  the 
composite  samples  gave  an  indication  of  the  validity  of  each  grab  sample 
as  the  run  progressed. 
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CHAPTER  V 

REVIEW  OF  PERTINENT  LITERATURE 

Mathematical  modeling  of  columnar  adsorption  reactions  for  design  and  opera- 
tion purposes  entails  the  description  of  the  mass  balance  for  the  system  and 
the  delineation  of  the  appropriate  equilibrium  and  kinetic  relationships. 

While  development  of  suitable  materials  balance,  equilibrium,  and  rate 
equations  is  simple  and  straightforward,  simultaneous  solution  of  this  set 
of  equations  is  mathematically  complex.  Frequently,  simplifying  assumptions 
and  numerical  techniques  are  required  to  enable  an  approximate  solution  of 
the  requisite  set  of  equations. 

Although  early  workers  (1)  showed  that  in  certain  cases  equilibrium 
data  alone  are  sufficient  for  such  mathematical  analyses,  it  it  generally 
necessary  to  consider  also  the  kinetics  of  the  sorption  process  to  permit 
accurate  calculation  of  the  breakthrough  curves.  Overall  rates  of  adsorption 
of  solutes  from  an  aqueous  stream  flowing  through  a fluidized  bed  of  granular 
active  carbon  represent  the  combined  effects  of:  (a)  diffusion  through 

the  hydrodynamic  boundary  layer  of  fluid  surrounding  the  porous  particle 
(film  diffusion);  (b)  diffusion  within  the  void  spaces  of  the  particle 
(pore  diffusion);  (c)  pore-surface  diffusion  along  the  walls  of  channels 
within  the  porous  sorbent  (solid-phase  diffusion);  and  (d)  phase-change 
in  the  pores  (adsorption).  Pore-surface  diffusion  is  distinguished  from 
pore  diffusion  in  that  it  occurs  after  phase-change  on  the  solid  side  of  the 
phase  boundary.  Under  different  operating  conditions,  any  combination  of 
these  resistances  may  be  operative.  Generally,  however,  only  one  or  two 
mechanisms  are  rate-limiting,  because  of  the  series  relationship  between  the 
various  resistances.  That  is,  the  slowest  rate  step  determines  the  overall 
rate  of  mass  transport.  For  sorption  on  granular  activated  carbon,  the  adsorp- 
tion step  itself  is  considered  to  be  extremely  rapid  (2).  Consequently,  only 
the  diffusional  resistances  are  of  import  in  determining  the  net  rate  of 


transport . 


Most  mathematical  solutions  for  systems  of  equations  describing  concentra- 
tion-time profiles  for  columnar  systems  are  limited  to  the  special  case  in 
which  only  one  of  these  four  rate  processes  controls  the  overall  rate  of 
adsorption.  Amundson  (3,4),  Edeskuty  and  Amundson  (5),  Hobson  and  Thcfos 
(6),  Hougen  and  Marshall  (7),  Kasten  and  Amundson  (8,9,10),  McCune  and 
Wilhelm  (11),  Resnick  and  White  (12),  Richardson  and  Szekely  (13),  Tien  and 
Thodos  (14),  Furnas  (15),  Goldstein  (16),  Klinkenberg  (17),  Kostecki , 

Manning  and  Canjar  (18),  and  Vermeulen  (19,20),  for  example,  have  presented 
solutions  for  systems  in  which  either  pore  diffusion,  film  diffusion, 
or  pore-surface  diffusion  is  the  only  significant  resistance  to  mass  trans- 
port. Several  investigators,  including  Colwell  and  Dranoff  (21),  Rosen 
(22,23),  Hall,  et  al , (24)  Masamune  and  Smi rh  1(24),  and  Tien  and  Thodos  (14) 
have  solved  the  more  general  case  in  which  both  trie  external  and  pore-diffusion 
or  pore-surface  diffusion  resistances  are  significant.  Masamune  and  Smith 
(25)  have  extended  these  concepts  and  presented  a solution  of  the  general 
three-resistance  case  in  integral  form.  Unfortunately,  however,  each  of 
the  above  developments  required  the  use  of  simplifying  assumptions  relative 
to  the  equilibrium  distribution  of  solute  between  solution  and  solid  phases. 

Of  the  four  basic  types  of  isotherms  --  irreversible,  favorable,  linear,  and 
unfavorable  --  each  respective  author  assumed  that  either  an  irreversible 
or  linear  distribution  represented  the  equilibrium  conditions,  primarily 
because  it  was  mathematically  expedient  to  do  so  to  enable  solution  of  the 
sets  of  partial  differential  equations.  Adsorption  isotherms  for  compounds 
typical  of  wastewater  pollutants  (26),  however,  are  generally  of  the 
favorable  type.  Certainly,  therefore,  equilibria  for  such  pollutants 
cannot  be  approximated  adequately  by  a linear  or  irreversible  type  of  isotherm. 
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If  approximations  to  the  kinetic  parameters  are  made,  such  as  the 
quadratic  driving-potential  approach,  solution  of  the  partial  differential 
equations  can  be  effected  for  non-linear  equilibrium  systems.  Hendricks 
(27),  Acrivos  (28),  Hall,  et  al , (24),  Lapidus  and  Rosen  (29),  and  Masamune 
and  Smith  (30,31),  used  various  approximations  for  the  rate  parameters  to 
obtain  solutions  for  breakthrough  profiles  for  the  case  of  favorable  equilibria 
of  the  Langmuir  type.  These  solutions,  however,  are  limited  to  the  condition 
when  only  one  of  the  possible  rate-limiting  mechanism  is  significant. 
Conversely,  the  second-order  reaction-kinetics  development  of  Thomas  (32,33) 
which  also  allows  the  use  of  a generalized  isotherm  of  the  Langmuir  type, 
provides  for  approximate  solutions  for  the  mono-,  d i - , and  tri-resistance 
cases.  Accordingly,  the  development  of  Thomas  appears  particularly  suited  for 
description  of  mass  transfer  processes  for  wastewater  systems.  Hiester  and 
Vermeulen  (34)  have  adapted  the  Thomas  solution  to  packed-bed  adsorbers,  and 
have  presented  methods  for  numerical  solution  of  the  equations  for  mass 
transfer  of  single-solutes  in  such  systems. 

Methods  by  which  rate  constants  and  mass  transfer  coefficients  have  been 
defined  by  those  investigators  who  provide  experimental  verification  of  their 
models  commonly  have  involved:  (a)  development  of  a model  based  on  several 
assumptions  regarding  equilibrium  and  rate  parameters;  (b)  assumption  of 
a succession  of  incremental  rate  constants  while  generating  a series  of  concen- 
tration-time profiles;  (c)  carrying  out  an  experimental  column  study  for  any 
suitable  system,  subject  to  the  nature  of  the  assumptions  involved;  and, 

(d)  determination  of  appropriate  rate  constants  for  the  system  by  means  of 
curve-matching  techniques.  It  is  to  be  noted  that  this  approach  is  not  truly 
predictive  in  nature,  as  it  must  resort  to  preliminary  experiments  of  a long- 
term pilot-plant  scale.  The  investigations  of  Allen,  et  al , (35,36),  provide 


an  excellent  example  of  one  such  approach  of  empirical  determination  of  rate 
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constants . 

One  of  the  principal  investigators  (37,38,39)  adapted  the  methods  of  Hiester 
and  Vermeulen  (34)  for  accurate  description  of  mass  transfer  processes  in 
columns  of  fluidized  active  carbon.  The  development  involved  correlation  of 
experimental  kinetic  and  equilibrium  data  for  agitated  non-flow  systems  with 
those  of  continuous  flow  fluid-bed  systems.  The  model  that  was  developed 
permits  prediction  of  concentration-time  profiles  for  systems  of  the  latter 
type  from  measurements  of  selected  adsorption  parameters  in  critically  designed 
non-flow  systems.  Measurements  in  simple,  non-flow  batch  experiments  can  of 
course  be  accomplished  much  more  readily  and  with  greater  ease  than  can  those 
in  relatively  complex  and  time  consuming  continuous  flow  column  tests.  Cookson 

(40)  recently  developed  a comparable  model  that  is  essentially  identical 
in  nature. 

While  the  studies  of  numerous  investigators  as  detailed  above  have 
provided  a model  which  serves  as  a rational  basis  for  the  design  of  fluidized- 
and  fixed-bed  adsorbers  for  the  removal  of  single  components  from  aqueous 
solution,  only  several  researchers  have  mounted  the  challenge  of  mathematically 
describing  adsorption  in  multicomponent,  competitive  systems  which,  of  course, 
are  most  commonly  encountered  in  prototype  applications.  Both  Collins  and  Chao 

(41)  and  Gariepy  and  Zweibel  (42)  developed  models  for  the  adsorption  of 
binary  mixtures  of  the  gas-phase  for  fixed-beds.  These  researchers  described 
competitive  adsorption  equilibria  either  by  the  competitive  Langmuir  model  or  by 
a thermodynamic , gas-phase  approach.  While  these  competitive  equilibrium 
expressions  were  relatively  satisfactory  for  carefully  selected  mixtures  of 
gases,  they  have  been  shown  to  be  unsatisfactory  for  description  of  competitive 
effects  for  adsorption  in  multicomponent  aqueous  solutions. 

Under  a previous  contract  the  principal  investigator  (43)  and  associates 
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developed  a mathematical  model  which  allows  simulation  of  the  performance 
of  adsorbers  when  treating  solutions  containing  multiple  components.  That 
development  employed  either  of  two  expressions  for  multiple-solute  adsorption 
equilibria;  that  of  Jain  and  Snoeyink  (44)  and  that  of  Radke  and  Prausni tz 
(45).  Moreover,  that  development  as  well  as  all  previous  models  that  describe 
the  performance  of  adsorption  from  multicomponent  solutions  in  columnar  contac- 
tors assumed  that  adsorption  on  granular  activated  carbon  was  fully  reversible. 
During  the  latter  portion  of  the  study  by  Keinath  (43)  this  assumption, 
however,  was  shown  not  to  be  entirely  valid.  Consequently,  there  remains  a 
demonstrated  need  for  modification  of  the  existing  model  to  account  for  all 
time-dependent  reversibi 1 i ty/irreversibi 1 i ty  considerations  such  that  mathematical 
simulations  of  adsorber  column  performance  accurately  predict  these  time  depen- 
dent phenomena.  Only  when  those  modifications  have  been  incorporated  into  the 
existing  model  will  it  be  generally  applicable  for  simulating  design  and  opera- 
tional alternatives  for  decontamination  of  multicomponent  aqueous  phase  systems 
by  use  of  adsorbers  which  contain  granular  activated  carbon  or,  for  that 
matter,  any  other  similar  adsorbent  material. 
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CHAPTER  VI 


MODEL  DEVELOPMENT  AND  IMPLEMENTATION 


MTERIAL  BALANCE  RELATIONSHIPS 


Because  adsorption  contactors  conventionally  are  of  a columnar  con 
figuration  and,  therefore,  have  concentration  profiles  in  the  axial 
direction,  it  is  necessary  to  conduct  a mass  balance  over  an  infinitesi 


mul  thickness  of  bed  at  a given  cross-section.  The  resultant  conserva 


solution  passing  through  that  section  must  equal  the  gain  of  solute  by 


the  adsorbent  contained  within  that  section.  For  the  purpose  of  this 


development  it  has  been  assumed  that  concentration  gradients  in  the  ra 


Verbally,  the  materials  balance  relationship  for  the  packed  bed 


case  may  be  expressed  as 


Input  to  Element  = Output  from  Element  + Adsorption  + Accumulation 


For  an  infinitesimal  thickness  of  bed  of  unit  cross-sectional  area  the 


following  mathematical  formulations  can  be  made  for  the  solution-phase 


Input  to  Element 


Output  from  Element 


Adsorption 


Accumulation 


U = solution  volumetric  flow  rate,  liters/hour 


C - solution-phase  concentration  of  solute,  moles/liter 


q = solid-phase  concentration  of  solute, moles/gram 
z = axial  distance,  cm 
t = time,  hours 

A = cross-sectional  area  of  column,  sq.  cm. 
l = void  ratio  or  porosity,  dimensionless 
i>  packed-bed  density,  grams/ liter 


Equating  in  accordance  with  Equation  VI-1  and  assuming  unit  cross-sec- 
tional area,  one  obtains: 


§dz; 


dz  + e 


‘DC 

at 


az 


which  upon  simplification  becomes: 


U — + 


0 


(VI-2) 


Solution  of  this  partial  differential  equation  (Equation  VI-2)  is  the 
essence  of  the  problem  of  mathematical  simulation  of  the  performance  of 
packed-bed  adsorption  contactors.  Such  solution  may  be  effected  by 
exact  solution  of  the  equations,  by  the  method  of  finite  differences, 
by  the  method  of  characteristics,  or  by  the  lumped  parameter  approach 
to  modeling  this  distributed  parameter  system.  To  obtain  an  exact  so- 
lution it  is  necessary  to  make  certain  simplifying  assumptions  regard- 
ing adsorption  equilibria;  such  as  a linear  or  irreversible  isotherm. 
Although  these  assumptions  are  relatively  restrictive,  they  did  permit 
the  development  of  exact  solutions  by  many  early  investigators  in- 
cluding notable  developments  of  Amundson  (4),  Kasten,  et  al . (10),  and 
Rosen  (23). 
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Solutions  by  the  methods  of  finite  differences  and  characteristics 
have  been  explicitly  detailed  by  Crank  (46)  and  Acrivos  (28),  respectively. 


Because  of  certain  complexities,  however,  these  are  not  user  oriented. 

In  contrast,  the  lumped  parameter  approach  is  user  oriented  as  it  permits 
the  use  of  the  continuous  system  modeling  program,  CSMPIII  , developed 
by  IBM  (47).  Consequently,  all  further  developments  relate  to  the 
lumped  parameter  approach. 


Lumped  Parameter  Approach 

Use  of  the  lumped  parameter  approach  requires  segmentation  of  the 


bed  of  adsorbent  into  a discrete  number  of  finite  elements  as  determined 


by  the  dispersion  characteristics  of  the  columnar  reactor.  Both  the 
liquid-  and  solid-phase  concentrations  of  solute  are  assumed  to  be  uni- 
form throughout  each  element.  Further,  continuity  of  mass  flows  of  the 
solute  between  adjoining  elements  must  be  maintained  for  both  the  li- 


quid and  solid  phases. 

For  the  case  in  which  an  adsorption  column  is  charged  with  an  ad- 
sorbent and  then  operated  in  the  unsteady-state  until  the  adsorbent  is 
entirely  exhausted,  the  following  formulations  can  be  made  for  an  ele- 
ment of  a packed-bed: 

Solution-Phase  Solid-Phase 


Input  to  Element: 
Output  from  Element: 
Adsorption: 
Accumulation: 


U * C 


U * C 


P * RA  * V 


dc  1 .. 
e n V 


«■ 
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where : 

V = volume  of  the  elements,  liters 
= rate  of  adsorption,  moles/gram-hour 
n » number  of  elements 

and  all  other  parameters  are  as  defined  above.  As  before,  equating  in 

accordance  with  Equation  VI-1  for  the  solution-phase: 

I dc  j 

U * C , U * C + r»  * K * V + r I — £ ' V 
n-1  n A |_  dt  ] 

which  yields  upon  rearrangement  for  the  n-th_  element, 
dc  u i ' 1 

JT  = rv|LCn-l  " Cn  I " lp/e  RA  (VI -3) 

Similarly,  for  the  solid-phase  the  mass  balance-relationship  for  the 
solute  becomes: 


0=0-p*Ra*V+p  (oqn/dt)  V 
which  simplifies  to, 

dqn/d<  = (VI-4) 

For  the  case  where  several  or  numerous  solutes  are  competing  for  the 
available  adsorption  sites,  mass  balances  for  the  solution-  and  solid- 
phases  must  be  written  for  each  solute  that  is  introduced  in  the  in- 
fluent to  the  column. 

Although  the  foregoing  equations  were  developed  for  packed-bed 
adsorbers,  they  can  readily  be  modified  for  the  fluidized-bed  mode  of 
operation.  This  involved  making  provision  for  mixing  of  the  adsorbent 
solids  within  the  column.  For  the  lumped  parameter  approach  this  is 

easily  accomplished  through  the  use  of  a procedural  solids  mixing  sub- 


It  is  important  to  recognize  that  the  set  of  differential  equations 
(Equations  VI-3  and  VI-4)  are  extremely  stiff.  That  is,  the  time  constant 
for  the  liquid-phase  equation  is  small  in  contrast  to  the  time  constant  for 
the  solid-phase  equation.  Numerical  integration  of  the  set  of  equations 
requires  specification  of  an  integration  interval  that  provides  for  the 
stable  and  accurate  integration  of  the  equation  that  responds  most  rapidly; 
i.e.,  the  liquid-phase  mass  balance. 

If,  however,  the  time  constants  for  the  two  equations  are  different  by 
several  orders  of  magnitude,  one  can  assume  that  the  equation  that  responds 
most  rapidly  is  continuously  at  steady-state.  For  this  case  the  term  dcn/dt 
in  Equation  VI-3  is  set  to  zero  and  the  resulting  algebraic  equation  is  solved 
for  Cn*  Specification  of  the  integration  interval  is  then  contingent  only 
on  the  remaining  differential  equation  - that  which  responds  most  slowly 
(solid-phase  mass  balance.  Equation  VI-4).  This  serves  to  materially 
decrease  the  computation  time  required  for  numerical  solution. 


ADSORPTION  EQUILIBRIA 


The  distribution  of  solute  between  liquid-  and  solid-phases  in  an 
adsorbent-solute-solvent  system  at  equilibrium  is  commonly  termed  an  ad 
sorption  isotherm.  Adsorption  equilibrium  data  is  conventionally  pre- 
sented and  correlated  by  plotting  the  quantity  of  solute  adsorbed  per 
unit  weight  of  adsorbent,  q,  as  a function  of  the  concentration  of  so- 
lute remaining  in  solution  at  equilibrium,  C. 

Several  mathematical  formulations  that  describe  adsorption  equili- 
bria and  that  have  been  widely  adopted  since  their  development  include 
the  model  isotherm  originally  proposed  by  Langmuir, 

_ Q + b * C 

q = TTT7T  (VI -5) 

that  of  Freundlich, 

q=K*C1/n'  (VI-6) 

and  that  of  Branauer,  Emmett  and  Teller, 


Q * A * C 

(Cs  - C)  (1  + LA  - lj  C/Cs) 


(VI-7) 


where : 

Q = ultimate  uptake  capacity  of  adsorbent,  moles/gram 
b = Langmuir  energy  term,  liters/mole 
K = adjustable  curve-fitting  constant 
n'=  adjustable  curve-fitting  constant 
A = BET  energy  term,  liters/mole 

Cs=  saturation  concentration  of  solute  in  solution,  moles/liter 


Of  these  three,  the  simple  empirical  Freundlich  expression  is  the  most 
widely  used  primarily  because  of  its  simplicity  and  because  it  has 
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been  successfully  employed  for  correlating  data  for  adsorption  of  solutes 
on  activated  carbon  over  limited  concentration  ranges.  The  Langmuir 
equation,  based  on  the  assumption  of  mono-layer  adsorption  on  a fixed 
number  of  equivalent  adsorptions  sites,  has  also  been  useful  over  limited 
concentration  ranges.  In  contrast,  the  BET  equation  which  is  based  on 
the  assumption  of  multiple-layer  adsorption  has  not  been  used  as  fre- 
quently as  either  the  Langmuir  or  Freundlich  models  principally  because 
of  the  difficulty  in  obtaining  suitable  values  for  Cs . 

LANGMUIR  COMPETITIVE  ADSORPTION  MODEL 

Each  of  the  three  models  which  describe  adsorption  equilibria, 
Equations  VI-5,  VI-6,  VI-7,  are,  of  course,  limited  for  application  single- 
solute systems.  Only  the  Langmuir  model  has  been  extended  to  account 
for  competitive  adsorption  equilibria  in  multiple-solute  systems.  For 
the  di-solute  case  the  extent  of  adsorption  for  solute  A is  given  by: 


n * b * r 
^A  ua 

qA  1 + bA  * CA  + bB  * CB 


(VI-8) 


and  for  solute  B the  solid-phase  equilibrium  concentration  may  be  ex- 
pressed as: 


= ^B  * bB  * CB 
^ 1 + bA  * CA  + bB  * CB 


(VI-9) 


where  the  constants,  Q^,  Qg,  b^,  and  bg,  are  those  that  are  measured  in 
mono-solute,  pure  solution  systems.  For  the  general  case,  where  (i)  is 
used  to  designate  the  number  of  solutes  in  the  system,  the  extent  of  ad- 
sorption for  solute  (j)  is  given  by: 


- W cj 

i 

1 + Z b,  * C, 
i=0 


(VI-10) 


Although  this  formulation  is  mathematically  simple,  it  has  been  experi- 
mentally verified  only  for  certain  selected  di-solute  competitive  sys- 
tems (44).  Accordingly,  it  is  important  to  recognize  that  this  relation- 
ship should  not  be  employed  indiscriminately  when  simulating  the  perfor- 
mance of  adsorption  contactors  for  competitive,  multiple-solute  systems. 
Rather,  the  relationship  must  be  experimentally  validated  for  the  speci- 
fic multi  pi  e-solute/sol  vent/adsorbent  system  before  its  use  can  be 
judged  appropriate. 


Semi-Competitive  Langmuir  Adsorption  Model 

Jain  and  Snoeyink  (44)  have  investigated  the  competitive  inter- 
action of  substituted  phenols  in  the  concentration  range  of  10"5  to 
10~2  moles  per  liter.  They  determined  that  the  competitive  Langmuir 
expression  was  not  satisfactory  for  description  of  results  obtained 
in  their  studies.  They  observed  a preferential  adsorption  for  one  com- 
pound in  comparison  to  the  other.  A fraction  of  the  material  which  was 
preferentially  adsorbed  appeared  to  respond  as  though  it  were  a single- 
solute system.  These  findings  led  to  the  development  of  the  following 
empirical  expressions 


\ ' V bA  CA  QbVa 

\ - \ * 6flCfl  1 + bACA  ♦ bBCB 

for  the  preferentially  adsorbed  material  and 

QBbBCB 

qB  1 + bACB  + bgCB 


(VI-11) 


(VI-12) 


for  the  less  readily  adsorbed  material. 
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The  first  term  of  Equation  VI-11  describes  that  fraction  of  the 

adsorbed  phase  which  is  preferent ial ly  adsorbed.  The  second  term 
describes  that  fraction  of  the  adsorbed  phase  that  is  competitive  in 
the  Langmuir  sense  with  the  second  solute.  The  rationale  for  the 
development  of  this  expression  was  based  upon  stcric  factors,  surface 
• ha  rue  I or i si  i c s , and  certain  physical  factors. 

I.qtiulion  II 1-1 2 represent  s the  competitive  interaction  of  tin- 
second  solute.  This  expression  is  exactly  the  same  as  that  of  the 
competitive  Langmuir  theory. 

Theoretically,  these  expressions  have  the  same  general  limitations 
as  the  Langmuir  equations.  Jain  and  Snoeyink  (44)  however,  have  shown 
that  these  equations  describe  a much  broader  spectra  of  competitive 
adsorption  isotherms  than  do  the  competitive  Langmuir  equations. 


Ideal  Solution  Theory  Adsorption  Isotherm 

Radke  and  Prausnitz  (45)  detailed  a procedure  for  the  prediction 
of  multi-solute  adsorption  froifi  dilute  solutions.  A detailed  descrip- 
tion of  their  method  is  given  below. 

The  thermodynamic  framework  of  the  Radke  and  Prausnitz (45)  method 
was  brsed  upon  the  solution-solid  interface  as  described  by  Gibbs. 
Assuming  the  solid  was  inert,  the  differential  equation  for  the 
Helmholtz  energy  of  the  adsorbed  phase  was  written  as 

dAa  = -SdT  - pdV  + Xdo  + V"uadna  + padna  (VI-13) 

4_  1 1 s S 


uhere  the  subscript  (i)  represents  solute  and  (s)  the  solvent.  The 
interfacial  tension  was  represented  by  X,  the  area  of  solution-solid 
interface  by  o,  and  moles  in  the  adsorbed  phase  by  na. 


Integration  of  liquation  VI-13  for  isothermal  conditions  and 
assuming  A,  pa , , V as  constant  gave 

Aa  = Aa  + gn>  + (A*  (VI-14) 

which  upon  differentiation  became  Gibbs  isothermal  adsorption  expression 

\do  ♦ ^n'‘d|d‘  a nad|,‘s'  - 0 . (VI-15) 

liit*  I i quid  - phase  which  was  not  affected  by  surface  forces  was 
expressed  as 


V'  1 , 1 1,1  n 

) n . du . + n . dp . =0 

/_■  t i li 


(VI-16) 


At  equilibrium  the  chemical  potential  of  the  liquid-phase  and  the  solid- 

cL  1 

phase  are  equal.  Accordingly,  multiplying  Equation  VI-1 6 by  ns/ns 


a , a _ 
+ n dp  =0 
s s 


(VI-17) 


n r \ . 

•y  Vnidu! 

n1  ^ 1 1 

s 

ami  subtracting  from  liquation  VI-15  the  expression  became 

adA  + Z!nidlii  ‘ "T  Zniddi  = ° ' (VI-18) 

ns 

Since  the  function  dA/dio  and  dA/dp^  are  the  Gibbs  excess  surface  func- 
tion, then  the  invariant  adsorption  (n?)  was  defined  as 

C . 

a 

- n.  - 

l 


m _ a l a 
n . - n . - a — n 
i C s 
s 


(VI-19) 


a a 

where  the  adsorbed  moles  of  solute  and  solvent  n‘.  and  n were  the  Gibbs 

1 s 

surface  excess  quantities.  Therefore,  the  following  expression 

-adA  = = odir  (VI -20) 

was  valid  over  the  entire  solute  range  where  the  spreading  pressure  (it) 
was  represented  by 
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" X (pure  solvent-solid)  - X (solution-solid)  (VI-21) 

Radke  and  Prausnitz  (45)  sliowcd  in  their  work  that  for  dilute 
solutions  n™  may  be  approximated  by 

nm  = VAC.  = q.  (VI-22) 

where  AC  was  the  decrease  in  concentration  or  constituent  in  a known 

i 

volume  ( V ) oi  the  solution-phase  when  contacted  with  a known  mass  of 
.ulsorheut.  It  was  Curt her  shown  that  it  was  useful  to  define  the  ad- 
sorbed phase  fugacity  as 

duf  = RT  d In  fa  (constant  T)  (VI-23) 

and 


lim 


fa 

1 


(VI-24) 


TT-^0  Z.7T 
1 

where  the  adsorbed  phase  mole  fraction  Z.  was  defined  as 
* i 


Z. 

i 


m 
n . 
i 

O m 

Lni 


m 

n. 

l 


n-r 


(VI-25) 


Considering  simultaneous  adsorption  of  the  solute  species  at  a 
constant  temperature  and  spreading  pressure,  it  was  assumed  that  an 
ideal  solution  phase  was  formed.  Based  upon  this  assumption,  the 
fugacity  f ^ must  be  proportional  to  the  mole  fraction  Z^ 

fi  = Z1  (VI-26) 

a° 

where  f‘.  is  the  fugacity  of  that  single-solute  in  dilute  solution  at 
the  same  temperature  and  spreading  pressure  of  the  mixture. 

liquation  VI-26  shows  that  n!J!  was  only  a function  of  the  adsorbed- 
phase  composition  Z.  and  invariant  adsorption  of  a single  solute  n°m, 


therefore,  n™  was  expressed  by 


HU  tfW 


into  liquation  VI-20  and  rearranging  the  terms 


ions 


which  upon  substitution  into  Equation  VI-28  yields  Equation  VI-27 


The  chemical  potential  of  solute  p' 


in  the  adsorbed-phase  must 


The  solid-phase  chemical  potential  was  expressed  as 


which  was  the  integrated  form  of  Equation  VI-23  . Since  T and  it 


established  the  single- solute  concentration  C 


following  substitution  for  p 


a was  made 


where  (*)  denotes  the  ideal-solution  standard  state  of  the  . 1 iquid-phase 
The  concentration  term  C?  (tt)  refers  to  a single-solute  adsorbed  from 


solution,  at  the  same  temperature  and  spreading  pressure  as  the  mixture 


ind  has  the  same  interpretation  as  n 


The  dilute  liquid  mixture  was  expressed  as 

u|(T,C)  = U-°(T)  + RT  In  C*  . (VI-33) 

Upon  substitution  of  Equations  VI-31  , VI-32  , and  VI-33  into  Equation 
VI-30  , the  following  expressions  were  obtained 

•-  f°(ji)Z^  (constant  T)  (VI-34) 

and 


C.[.Xi  = C°(tt)2.  (constant  T)  . (VI-35) 

lhe  spreading  pressure  must  be  evaluated  before  the  mixture  con- 
centration and  adsorbed  phase  may  be  computed.  The  spreading  pressure 
may  be  computed  by  integration  Equation  VI-20  to  give 


C?  m° (r  . 
1 n.  (C.) 


T /'in.  u 

ir(C?)  = — / — — 

1 0 J C° 


dC° 

l 


(VI-36) 


0 i 

Experimental  single-solute  isotherms  are  used  to  develop  a curve  of 
m°  o o 

n.  /C.  as  a function  of  C. . These  can  be  either  numerically  or 
ii  l 

graphically  integrated  to  determine  n.  A graphical  procedure  of  Myers 
and  Prausnitz  (48)  and  a numerical  method  have  been  described  in  detail 
in  Appendix  B.  Because  the  single-solute  data  describes  7r,  there  is 
no  need  for  a theoretical  model  of  simplifying  or  restrictive  assump- 
tions to  describe  the  adsorption  equilibria. 

Multi-solute  isothermal  equilibria  may  be  calculated  using  a 
method  described  by  Myers  and  Prausnitz  (48)  . The  procedure  is  sum- 
marized below: 

(i)  Experimental  single-solute  isotherm  data  is  integrated 
according  to  liquation  VI-36  to  yield  the  spreading 
pressure  for  each  solute  as  a function  of  concentration. 


This  is  given  by: 
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"j  = fjCC?)  , and 


(VI-37) 


"2  = MCP  • 


(ii)  The  equilibrium  relationships  are  given  by: 

C.j.Xj  = C°Zj  , and 

c,.x2  = C-Cl  - Zl)  . 


(VI-38) 


(VI-39) 

(VI-40) 


Those  are  shown  in  Figure  VI-1. 

(iii)  The  spreading  pressure  of  each  constituent  is  equivalent 
to  that  of  the  mixture,  i.e. 


IT,  = TT  = TT  . 
1 2 


(VI-41) 


(iv)  For  given  values  of  and  xj  the  remaining  variables 

it  , at 2 , z C° , and  C°  may  be  determined  from  the  preceding 
five  equations. 

(v)  Using  the  single-solute  data  and  Equation  VI-27,  the  total 
invariant  adsorption  may  be  evaluated.  The  invariant  ad- 
sorption of  each  species  in  the  mixture  may  be  determined 


by  the  following  balances 


m m,  , 

ni  = n.j.Z  , and 


m rn  f . -7  > 

n2  = nT(l  - Z:)  . 


(VI-42) 


(VI-43) 


The  ideal  solution  theory  provides  a method  of  predicting  the 
adsorption  of  two  or  more  solutes  from  solution.  This  method  is  based 
upon  sound  thermodynamic  relationships  and  there  are  no  restrictions  on 
miscibility  nor  mono-layer  theory.  Considering  the  restrictions  of  the 
competitive  Langmuir  as  well  as  semi-competitive  Langmuir  developments, 
the  ideal  solution  theory  method  has  the  most  fundamental  basis- 


i 
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KINETICS  CF  INTERPHASE  SOLUTE  TRANSPORT 

To  relate  the  dynamically  varying  solid-  and  solution-phase  concen- 
trations of  solute  at  any  time  and  position  in  an  adsorption  contactor, 
it  is  necessary  to  fully  understand  the  kinetics  of  solute  transport. 
Consequently,  the  mechanisms  of  such  transport  must  be  considered  in 
detai 1 . 

The  overall  rate  of  adsorption  of  solute,  from  a solvent  stream 
flowing  through  a bed  of  porous,  granular  adsorbent  represents  the  com- 
bined effects  of:  (1)  diffusion  through  the  boundary  layer  of  fluid 

surrounding  the  adsorbent  particle  (film  or  external  diffusion);  (2) 
diffusion  within  the  pores  of  the  particle  (pore  diffusion);  (3)  diffu- 
sion along  the  surface  of  the  pores  (pore-surface  diffusion);  and  (4) 
adsorption  on  the  internal  pore  surfaces  (adsorption).  Pore-surface 
diffusion  is  distinguished  from  pore  diffusion  in  that  it  occurs  on 
the  solid  side  of  the  phase  boundary.  That  is,  the  solute  may  enter  a 
particle  of  adsorbent  directly  from  the  exterior  surface  by  movement  in 
a condensed  or  adsorbed  layer  along  the  pore  surfaces,  or  it  may  diffuse 
through  the  fluid-phase  held  in  the  pores  and  then  be  deposited  at  a 
stationary  location  on  the  pore  surface.  In  pore  diffusion,  the  mass 
transfer  process  occurs  before  the  phase  change,  while  in  pore-surface 
diffusion  it  occurs  afterward.  Under  different  operating  conditions, 
the  same  particle  can  exhibit  either  surface-  or  pore-diffusion  behavior, 
favoring  the  former  if  the  concentration  level  in  the  pore-fluid  is  low 


sion/pore-surface  diffusion/adsorption,  the  slowest  of  the  three  mass 
transfer  mechanisms  will  in  each  case  limit  the  overall  rate  of  adsorp- 
tion because  these  rate  steps  occur  in  a series  relationship. 
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Adsorption  Rate  Step 

It  has  been  experimentally  determined  that  surface  tension  attains 
equilibrium  after  a disturbance  in  approximately  ten  milliseconds. 

This  circumstantial  evidence  indicates  that  the  adsorption  process  it- 
self is  a relatively  rapid  process  and,  therefore,  is  probably  not  rate 
controlling.  Considerable  experimental  evidence  exists  that  supports 
this  supposition.  Accordingly,  adsorption  is  generally  considered  as 
not  limiting  the  overall  rate  of  adsorption.  It  will,  therefore,  not 
be  considered  further  in  this  development. 


Pore-Diffusion 

For  a porous  adsorbent  the  rate  of  diffusion  into  the  liquid  in 
the  pores  of  a spherical  particle  is  given  by: 


D 


pore 


~52c  2 9C 

3r2  r H 


3t 


(VI-44) 


where: 

2 

„ = pore  diffusivity,  cm  /sec 
pore 

r = radial  distance  within  particle,  cm 

X = internal  porosity  of  adsorbent  particles,  dimensionless 
Pp  = density  of  adsorbent  particle,  gm/liter 

The  mean  concentration  of  the  entire  particle  of  diameter,  Dp  (cm),  is 
obtained  by: 


0J2 

r 


q r dr 


(VI-45) 


| 


I 
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For  liquids,  the  pore  diffusivity  term  (Dpore)  may  be  approximated  by 
the  relationship: 


D = Di  *' 
pore  — s — 


(VI-46) 


where: 


D = diffusivity  of  solute  in  solvent,  cm  /sec 


Although  Equations  VI-44,  VI-45,  and  VI46  clearly  describe  the  pore- 
diffusion  rate-limiting  case,  it  is  to  be  noted  that  solution  of  these 
equations  in  conjunction  with  either  Equations  VI-2  or  VI-3  (materials 
balance  relationships  for  the  columnar  adsorption  contactor)  has  not  yet 
been  accomplished,  either  exactly  or  numerically.  This,  of  course,  is 
due  to  the  fact  that  an  additional  independent  variable,  radial  distance 
within  the  porous  adsorbent  particle,  is  introduced  into  the  mass  balance 
relationship  for  a packed-bed  contactor.  Consequently,  a linear  driving 
potential  is  frequently  employed  as  a gross  approximation  (49).  The  rate 
equation  then  is: 


R 


= = 
A dt 


k * A * (C  - C*) 


(VI-47) 


where: 


kp  = mass-transfer  coefficient  (pore),  cm/hr 

? 

Ap  = external  interfacial  transfer  area  for  adsorbent,  cm  /liter 
C*  = solution-phase  concentration  of  solute  considered  to  be  in 
equilibrium  with  the  outer  surface  of  the  adsorbent  particle 
moles/liter. 

For  this  approximation,  the  term  kp  * Ap  is  evaluated  by: 


60  * D „ 

k * a = P0|,'e  (1  _ r\ 

> P D 2 U £) 


(VI-48) 
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Pore-Surf aoe  Diffusion 


Diffusion  of  solute  molecules  along  the  surfaces  of  the  pores  in 
a porous  adsorbent  is  described  by: 


n i!a+  lia  = 
s ar2  r 3r  3t 


where: 

2 

Ds  = pore-surface  diffusivity,  cm  /sec. 

As  before,  solution  of  Equation  VI-49  in  association  with  Equation  VI-2 
or  VI-3  cannot  currently  be  achieved  with  the  various  numerical  methods 
that  are  available.  This  equation,  therefore,  is  frequently  approximated 
by  the  1 inear-dri ving  force  relation  of  Glueckauf  and  Coates  (50): 

RA  = = ks  * Ap  (q*  ‘ q)  (VI'50) 

where: 


= mass-transfer  coefficient  (pore-surface) , cm/hr 

q*  = solid-phase  concentration  of  solute  considered  to  be  in  equili- 
brium with  the  instantaneous  fluid-phase  concentration  outside 
the  particle,  moles/gram 

The  product  * Ap  is  related  to  the  diffusivity  and  to  the  adsorbent 

particle  diameter  through  the  equation: 

60  * D. 


ks  * Ap  = 


D. 


(VI-51) 


Experimentally  determined  values  for  D$  are  required  for  solution  of 
Equation  VI-51  such  that  estimates  of  the  term  Kg  * Ap  can  be  obtained. 
Methods  of  obtaining  values  for  D$  have  been  detailed  explicitly  by 
DiGiano  (51).  For  aqueous  solutions  of  various  substituted  phenols 
and  benzene  sulfonates,  DiGiano  (51)  determined  ratios  of  D /Df  which 
varied  from  0.1  to  0.3. 
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A somewhat  better  approximation  to  the  behavior  of  Equation  VI-49  is 
given  by  the  quadratic  driving  potential  form  postulated  by  Vermeulen 

(19): 


= = 
dt 


jlJL 


2q 


(VI-52) 


where : 

= initial  solid-phase  concentration  (usually  zero),  moles/gram 
v = 1 / ( R + 15  (1  - R)/tt2) 

R = separation  factor,  dimensionless 

Film  or  External  Diffusion 

Transfer  of  solute  molecules  from  the  bulk  solution  through  the 
hydrodynamic  boundary  layer  surrounding  the  granular  adsorbent  particle 
to  its  external  surface  is  commonly  described  by: 

",'3?  ’ kf  * \ <c  - c*>  (VI-53) 

where: 

k^  = film-diffusion  controlled  mass  transfer  coefficient,  cm/hr 

Use  of  Equation  VI-53  does  not  necessarily  imply  the  existence  of  a stag- 
nant hydrodynamic  boundary  layer  through  which  molecular  diffusion 
occurs.  Rather,  this  equation  may  generally  be  used  for  description 
of  mass  transfer  through  the  boundary  layer.  Certainly,  this  would 
include  the  effects  of  eddy  transport  in  addition  to  molecular  diffu- 
sion. 

Although  various  theories  including  the  penetration  and  free  sur- 
face models  have  been  advanced  to  describe  transient  diffusion  through 
an  external  boundary  layer,  the  straightforward  relationship  given  by 

Equation  VI-53  is  the  only  form  that  is  generally  applicable  for  simula- 
tion of  the  dynamics  of  a columnar  adsorber.  Because  of  the  difficulty 


in  describing  the  hydrodynamic  conditions  in  the  boundary  layer,  mass 
transfer  coefficients,  k^  are  not  determined  directly  but  are  obtained 
by  correlation  of  the  mass  transfer  factor,  j^,  with  the  Sherwood  num- 
ber, Reynolds  number,  and  Schmidt  number  as  follows: 


Jc!  ■ NSh/(N 


, * n 

’Re  NSc  ‘ 


(VI-54) 


where: 


= (k^/U)  * (v/D„)^'/^,  mass  transfer  factor,  dimensionless 

NPU  = k,  * D /D„,  Sherwood  number,  dimensionless 
Sh  f p 2 

^Sc  = » Schmidt  number,  dimensionless 

NRe  = U * D /v,  Reynolds  number,  dimensionless 
U = U/A,  velocity  of  flow,  cm/hr 

2 

v = kinematic  viscosity  of  solution,  cm  /sec 


Numerous  investigators  have  conducted  intensive  experimental  studies 
to  evaluate  the  correlation  given  in  the  foregoing  equation.  Equation 

VI-54,  for  a variety  of  different  system  types  includinq  qas-liquid,  li- 
quid-solid, and  gas-solid  two-phase  systems.  Extensive  investigations 
have  also  focused  on  single-particle  as  well  as  multiparticle  systems 
which  are  operated  either  in  the  packed,  expanded,  semi fluidized,  or 
fully  f luidized-bed  modes  (52,  53,  54,  & 55). 

Accordingly,  if  it  has  been  determined  that  film  diffusion  or  ex- 
ternal transport  is  rate  limiting,  then  one  can  obtain  a suitable  mass 
transfer  coefficient  from  one  of  the  existing  correlations.  It  is  im- 
portant, however,  to  pay  particular  attention  to  the  type  of  system  for 
which  the  correlation  was  obtained  and  to  the  procedure  employed  for 
calculation  of  mass  transfer  coefficients. 


MODEL  IMPLEMENTATION 


Prediction  of  the  performance  of  adsorption  contactors  involves 
the  solution  of  a set  of  differential  equations  which  describe  the  fluid- 
and  solid-phase  material  balances.  The  computer  provides  a convenient 
method  of  structuring  the  differential  equations  for  numerical  solution. 
Structure  of  the  model  may  be  directed  toward  the  use  of  any  number  of 
numerical  integration  procedures.  The  integration  technique  selected 
depends  upon  the  users  familiarity  with  the  various  techniques  as  well 
as  convergence  criteria  of  the  technique  selected.  For  this  reason  it 
was  decided  to  structure  the  model  such  that  the  user  might  select  the 
integration  routine  directly  from  the  scientific  subroutines  stored  in 
the  computer  library.  This  would  permit  the  engineer  who  has  minimum 
programming  experience  to  utilize  and  understand  the  function  of  the  model, 
while  allowing  the  engineer  who  is  more  adept  in  computer  programming 
the  versatility  of  restructuring  the  model  to  any  other  numerical  technique 
he  desired. 

For  this  reason  the  IBM  System/360  Continuous  System  Modeling  Program 
(CSMPIII)  was  used  for  solution  of  the  differential  contacting  model.  This 
system  is  a problem  oriented  program  which  permits  the  problems  to  be 
prepared  directly  from  block  diagrams  or  a set  of  ordinary  differential 
equations.  It  contains  a basic  set  of  34  functional  blocks  which  may  be 
used  to  represent  components  of  a continuous  system.  This  basic  set 
of  functional  blocks  includes  conventional  analog  computer  components  such 
as  integrators,  relays,  and  time  delays  as  well  as  specialty  functions  such 
as  logic  functions  and  macro  systems.  In  effect  CSMP  III  does  not  have  to 
operate  within  the  simulation  language  framework  but  may  be  structured 
in  the  symbols  characteristic  of  the  language  of  any  specialty  field. 

The  program  will  also  accept  FORTRAN  statements,  thereby  allowing  the  user 
to  solve  complex  nonlinear  and  time-variant  problems. 


r* 
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As  was  suggested  above,  solution  of  the  set  of  differential  equations 
that  describe  the  mass  balances  for  the  solution-  and  solid-phases  of  an 
adsorption  contactor  is  facilitated  by  use  of  the  distributed  parameter 
approach  when  used  in  conjunction  with  CSKP  III.  Programs  formulated  in 
CSMP  III  generally  consist  of  three  segments:  INITIAL , ’DYNAMIC , and 

TERMINAL.  Broadly,  the  INITIAL  segment  is  intended  exclusively  for 
computation  of  initial  condition  values.  The  DYNAMIC  segment  is  normally 
the  most  extensive  of  the  model.  It  includes  the  complete  description  of  the 
system  dynamics,  together  with  any  other  computations  desired  during  the 
run.  Conversely,  computations  which  are  to  be  performed  only  after  comple- 
tion of  the  dynamic  run  are  placed  into  the  TERMINAL  segment  of  the  model. 
These  will  often  be  simple  calculations  or  optimization  algorithms  based 
on  the  final  value  of  one  or  more  model  variables. 

For  an  adsorption  contactor  the  differential  equations,  VI-3  and  VI-4 
describe  the  system  dynamics  or  the  time-variant  performance  of  each 
segment  of  the  adsorbent  bed.  Of  course,  these  require  input  of  the  appropri- 
ate rate  relationship  for  R^  (Equations  VI-47,  VI-50,  VI-52,  or  VI-53) 
and  a suitable  description  cf  the  competitive  equilbrium  distribution  of 
solute  between  solution  and  solid  phases.  If  it  is  assumed  that  mass-transfer 
is  rate  limiting  and  that  the  competitive  Langmuir  equation  adequately  re- 
presents the  pertinent  adsorption  equilibria,  then  the  system  dynamics  can 
be  formulated  as  shown  below  in  the  DYNAMIC  segment  of  CSMP  III. 


DYNAMIC 


♦**  DIFFERENTIAL  EQUATIONS  FOR  ADSORPTION  BED 


*** 


QA* INTGRL ( I CQA*  QADOT • 1 1 I 
Q6*INTGRL(I CQS*  08 OD  T * 1 1 I 


V.  . — - « _ ■ _ _ 
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PROCEDURE  QADOT, QB30T*ADS I QA, OB, I SOB, QAMAX * BA, KFA, KFB, A I 

oa  10  1=2,11 

IF  I OBI  I ) .LT.5.E-5)  OB  t MV=2.E*A 
IF  ( QB( I I . GE . 5 . E-5 ) QBl NV  = 1/OBI  I ) 

C8INV=FUNGEN( I SOB, 5, OBI  NV) 

8SL0PE  = SL3PE II  SOB, 3, OBI  MV) 

QBMAX=1/ (QBINV-I l / BSLOPE ) *C B I NV ) 

88=8SL0PE/QBMAX 

C$TARA=QRMAX*BB*QA( I ) / I OAM AX*QB* AX*BA*BB -QAMAX *QB II ) *BA*B8, . . 
-QBMAX*QAI I ) *8A*BB) 

CSTARB=QAMAX*BA*QB(  I )/<QAMAX*QB»AX*8A*BB-QAMAX*QB( I )*BA*88... 
-OBMAX*OAI I ) *BA*8B ) 

C A I I ) = I DUHl *CAI  I-l)*DUM2*CSrARA)/IDUMUDU«2» 

CB ( I )=(0UMl*C8<  l-l)«-DUM3*CSTARB)/(OUMl*OUM3) 

3A03TI I >=KFA*A*ICAI I ) -CSTAR A ) • l . E-3 
OBDOT I I )=<FB*A*I CB1 I 1-CSTARB )*l.E-3 
10  CONTINUE 

ENDPROCEDURE 

Initial  conditions  for  second  and  fourth  equations  of  each  element 
of  the  bed  are  generally  specified  in  the  INITIAL  segment  of  the  CSKP 
III  program  using  an  INCCN  or  TABLE  label  as  shown  in  the  program  listing 
which  follows: 

***  INITIAL  CONDITIONS  FOR  SOLID  PHASE  CONCENTRATION  *** 

************************************************************************ 


TABLE  IC0AI2-11 1*11*0.0 
TABLE  ICQBI2-ll)*ll*0.0 

Likewise,  parameters  which  must  be  specified  for  the  simulation  are 
introduced  in  the  INTIAL  segment  by  use  of  the  PARAN:  statement  as  may 
be  observed  in  the  listings  below: 


I 
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PARAMETER  INPUTS  REQUIRED  ARE  DEFINED  AS  FOLLOWS  *** 

********************************************************************* 


#** 

CARBON  MASS  OF  CARBON  IN  BED  (GRAMS)  *** 

A EXTERNAL  TRANSFER  AREA  FOR  CARBON  (SQ  CM/GRAM)  *** 

DP  DIAMETER  OF  CARBON  PARTICLES  (CM)  *** 

RHOPAC  BACKED  3ED  DENSITY  (GRAMS/LITER)  *** 

E PS PACKED  BED  POROSITY  (DIMENSIONLESS) 

QAMAX  LANGMUIR  ULTIMATE  UPTAKE  CAPACITY  FOR  A (MOLES/GRAM)*** 

QBMAX  LANGMUIR  ULTIMATE  UPTAKE  CAPACITY  FOR  B ( MOL ES/ GRAM ) *** 

BA  LANGMUIR  ENERGY  TERM  FOR  A (LITERS/MOLE)  *♦* 

B8  LANGMUIR  ENERGY  TERM  FOR  B (LITERS/MOLE)  *** 

C A ( 1 ) SOLUTION-PHASE  CONCENTRATION  OF  A ( MOLES/L I TER  I *** 

CB(l)  SOLUTION-PHASE  CONCENTRATION  OF  B (MOLES/LITER)  **♦ 

OLA  DIFFUSIVITY  OF  SOLUTE  A (SQ  CM/SEC)  *** 

DLB  DIFFUSIVITY  OF  SOLUTE  B ( SQ  CM/SEC)  *** 

AREA  CROSS-SECTIONAL  AREA  OF  COLUMN  (SQ  CM)  *** 

U SOLUTION  VOLUMETRIC  FLOW  RATE  (LITERS/HOUR)  *** 

GAMMA  KINEMATIC  VISCOSITY  OF  WATER  (SQ  CM/SEC)  *** 

EXP RATIOS  EXPANDED  /PACKED  BED  VOLUME  ( OIMENS IONLESS  ) *** 

*** 


********************************************************************* 
PARAMETER  INPUTS  FOLLOW  *♦* 

********************************************************************* 

PARAM  CARB0N=36.79 
PARAM  A*12*.6 
PARAM  DP=9.55E-2 
PARAM  RHOPAC=253. 

PARAM  EPS=0.<f3 
PARAM  QAMAX=3.4958E-3 
PARAM  BA=653 . 

TABLE  C A ( l ) * l . E-3 
TABLE  CBII)*L.E-3 
PARAM  CA0=l.£-3 
PARAM  CBO=l.E-3 
PARAM  DL A=8. 672E-6 
PARAM  DLB=8.5AE-6 
PARAM  AREA**. 978 
PARAM  U-6.12 
PARAM  GAMMA*9.3A8E-3 
PARAM  EXP* ( l • 0* l • 5 ) 


Moreover,  certain  pararreters  which  are  non-time-variant  and  must  be 
calculated  prior  to  the  simulation  are  evaluated  in  the  INITIAL  segment. 

A listing  of  those  parameters  that  rr.u s t be  determined  by  calculation  follows: 


************************************************* ************** ********* 

***  CALCULATED  INPUT  PARAMETERS  ARE  DEFINED  AS  FOLLOWS  *** 

************************************************************************ 


***  *** 

***  RHO EXPANDED  BED  DENSITY  I GRAMS/  LITER)  *** 

***  VP  TOTAL  VOLUME  OF  PACKED  BED  (LITERS)  *** 

***  V TOTAL  VOLUME  OF  PACKED  BED  (LITERS)  ♦♦♦ 

***  EPSI  EXPANDED  BED  POROSI TY  (DIMENSIONLESS)  *** 

***  UBAR  SUPERFICIAL  VELOCITY  OF  FLOW  THROUGH  BED  (CM/SEC)  *** 

***  NREMOD  MODIFIED  REYNOLDS  NUMBER  (DIMENSIONLESS)  *** 

***  NSC  A SCHMIDT  NUMBER  FOR  A (DIMENSIONLESS)  *** 

***  NSCB  SCHMIDT  NUMBER  FDR  B (DIMENSIONLESS)  *** 

M*  JO MASS  TRANSFER  FACTOR  (DIMENSIONLESS)  **• 

***  KFA  MASS  TRANSFER  COEFFICIENT  FOR  A (CM/HOUR)  *** 

***  KFB  MASS  TRANSFER  COEFFICIENT  FOR  B (CM/HOUR)  *** 

***  DUM1  DUMMY  COEFFICIENT  *** 

***  DUM2  DUMMY  COEFFICIENT  *** 

***  DUM3  DUMMY  COEFFICIENT  *** 

***  *** 


**** ******* ****************** *** **************************************** 
***  INPUT  PARAMETER  CALCULATIONS  FOLLOW  *** 

************************************************************************ 


RHO=RHOPAC/EXP 
VP3CARBON/{ R HOP AC *10. ) 

V*C ARBON/ ( RHO*IO.  ) 

EPSI  = (EPS*VP*-V-VP)/V 

UBAR=(U*0. 277777) /AREA 

NREMOD=(UBAR*DP)/(GAMMA*( I. -EPS!  ) ) 

NSCA=GAMMA/DLA 

NSCB= GAMMA /DLB 

JD=l.34*NREM30**( -0.468) 

KFA=< JD*JBAR/NSCA**( 2./ 3. ) )*3600. 
KFB=(JD*UBAR/NSCB**(2./3.  ) )*3600. 
DUM1=*U/(EPSI*V) 
DUM2»(KFA*A*RH0*l.E-3)/EPSI 
DUM3»(KF8*A*RHO*l.E-3)/EPSI 


If,  in  constrast,  certain  parameters  are  time  dependent,  they  must  be  included 
in  the  DYNAMIC  segment  of  the  model,  such  that  at  any  time  during  the  integra- 
tion the  parameter  is  evaluated  for  its  current  numerical  value. 

In  addition  to  the  foregoing  listing  in  which  the  dynamic  model  is  fully 
specified,  it  is  necessary  also  to  provide  for  several  execution  and  output 
control  statements.  Execution  control  statements  are  used  to  specify  certain 
items  relating  to  the  actual  simulation;  e.g.  integration  method,  run  time,  in- 
tegration interval,  relative  error,  output  times,  etc.  An  example  of  the  execu- 
tion control  statements  used  in  the  current  simulation  studies  is  shewn  below. 
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************************************************************************ 
***  EXECUTION  CONTROL  STATEMENTS  *** 

*************** **************************************************  ******* 


METHOD  RECT 

TIMER  DELT=.?n'F I NT I M = 48. 0 , PRDE L = 0. 2 , 0UTDEL=0.  A 


The  first  of  two  control  statement  cards  specifies  the  type  of  integration 
method  which  is  to  be  employed  during  the  simulation.  For  this  case  the 
Rectangular  Rule  (RECT)  integration  method  was  utilized.  Other  integration 
methods  which  are  available  are  listed  in  the  CSMP  III  User's  Manual  (47). 

On  the  Timer  statement  one  specifies  the  integration  interval  (DELT),  the 
total  duration  of  the  simulation  run  (FI NT I M ) , and  simulation  data  output 
intervals  (PRDEL  & OUTDEL). 

For  the  simulations  conducted  as  part  of  this  study  the  output  was 
obtained  in  a graphical  form  in  the  format  presented  by  either  a line-printer 
or  x-y  plotter.  An  example  of  the  output  control  statements  used  in  the 
current  simulation  studies  is  shown  below: 

************************************************************************ 
***  OUTPUT  CONTROL  STATEMENTS  *** 

************************************************************************ 

PRINT  QA(2-lll .QBI2-III 
OUTPUT  TIME«CA(ll),CB(lU 
PAGE  GROUP* (0.,2.E-3) 

OUTPUT  TIME, OAT, QBT 
PAGE  GROUP =( 0.,2.E-3I 
OUTPUT  TIME,CA(li),CB(ll) 

PAGE  XYPLQT»HEIGHT=5.0»WIDTH=8.0» GROUP* (0«,2.E-3I 
OUTPUT  TIME, OAT, QBT 

PAGE  XYPL3T.HE IGHT*5.0,WIDTR*B.O,OROUP*<0.,2.E-3) 

********** ******************************************************* ******* 

To  account  for  the  mixing  of  the  granular  activated  carbon  in  fluidized 

beds  a procedural  solids  mixing  subrouting  is  required.  The  subroutine  employed 
for  the  current  studies  is  listed  as  follows: 


*♦******♦*********************♦*♦*♦*♦♦**♦***♦*♦♦♦♦*♦♦•♦♦*♦**♦♦****♦**♦*♦ 

***  TO  ACCOUNT  FOR  MIXING  OF  CARBON  IN  FLUIDIZED  BEDS  (EXPM.O),  *** 

***  AVERAGE  SOL  I D-PHASE  CONCENTRATIONS  ARE  CALCULATED  AND  QA'S  L **• 

***  OB'S  ARE  REINITIALIZED  AT  THE  AVERAGE  VALUES  *** 

************************************************************************ 

PROCEDURE  3AT,QBT=SOR3EDIQA,QB) 

QATDT  = 0 . 

QBTOT  *0 . 

DO  1L  J = 2 » 1 1 
QAT3T*0AT0T*3A( JI 

11  QBTQT=QBT3T*QB(J ) 

QAT  = QATOT/lrJ. 

QBT=QBT0T/1 0. 

IF  (EXP. E3. 1.0)  GO  TO  13 
DO  12  K = 2 » 1 1 
QA ( K ) =QAT 

12  QB(K)=QBT 

13  CONTINUE 
ENDPPOCEDURE 

************************************************************************ 


A complete  listing  of  a typical  program  employed  for  the  simulations 


has  been  attached  as  Appendix  C. 
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CHAPTER  VII 
EXPERIMENTAL  RESULTS 

ADSORPTION  EQUILIBRIA  STUDIES 

To  obtain  equilibrium  parameters  for  use  in  the  basic  mathematical 
model  for  simulation  of  effluent-time  profiles  for  columnar  flow-through 
systems,  it  was  necessary  to  obtain  single-solute  adsorption  isotherms. 

These  were  obtained  for  three  organic  herbicides;  dinitro-o-sec-butylphenol 
(DNOSBP),  o-phenyl phenol  (OPP),  and  2 ,4-di chlorophenol  (2,4-DCP).  The 
isotherms  are  shown  in  Figures  V 1 1 - 1 , VII-2,  and  V 1 1 - 3 , respectively. 

ANALYSIS  OF  COMPETITIVE  ADSORPTION  EQUILIBRIA 

All  columnar,  continuous-flow,  bisolute  competitive  adsorption  experi- 
mental data  was  analyzed  with  respect  to  the  steady-state  condition  (equilibrium) 
which  occurred  at  the  end  of  each  experimental  run.  The  steady-state  values 
were  compared  with  multi-solute  adsorption  equilibria  predictions  of  the 
Langmuir  competitive  and  semi -competitive  models  (previously  described) 
and  with  a graphical  approach  termed  the  Graphical  Irreversible  Model. 

The  graphical  model  is  analogous  to  the  Semi -Competiti ve  Langmuir 
Model  in  that  it  accounts  for  a preferentially  adsorbed  solute  and  a competi- 
tively adsorbed  solute  -in  bisolute  systems.  The  method  is  based  on  the 
assumption  that  the  maximum  quantity  of  both  solutes  that  can  be  adsorbed  is 
equal  to  that  quantity  of  the  preferentially  adsorbed  solute  that  would  be 
adsorbed  if  it  were  applied  as  a single-solute  at  the  reference  solution- 
phase  concentration. 

The  Graphical  Irreversible  Model  is,  as  might  be  expected,  based  on  the 
assumption  of  irreversibility;  that  is,  once  a solute  is  adsorbed,  it  remains 
adsorbed  and  cannot  be  displaced  by  another  solute.  The  method  can  best  be 
illustrated  by  example.  If  an  adsorber  were  initially  brought  to  equilibrium 
with  DNOSBP  at  a solution-phase  concentration  of  80  mircomoles  per  liter. 


FIGURE  VII-1 : ADSORPTION  ISOTHERM  FOR  DNOSEP 
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then  the  corresponding  solid-phase  concentration  should  attain  a value  of 
700  micromoles  per  gram  (Figure  VI I -1 ) . If  OPP  were  subsequently  introduced 
to  the  adsorber  along  with  the  DNOSBP  at  a concentration  of  40  micromolar,  then 
the  solid-phase  concentration  of  OPP  should  attain  a value  of  approximately 
200  micromoles  per  gram.  This  is  determined  by  establishing  the  quantity  of 
OPP  that  would  be  absorbed  if  it  were  introduced  to  the  adsorber  as  a single- 
solute at  the  40  micromolar  level  (ca.  900  micromoles  per  gram  as  noted  in 
Figure  VII-2)  and  then  subtracting  that  quantity  of  the  competitive  solute 
that  had  already  been  adsorbed  (700  micromoles  per  gram). 

Conversely,  if  OPP  had  been  applied  to  the  adsorber  as  a single-solute 
at  the  40  micromolar  concentration  level  before  DNOSBP  was  added  to  the 
influent  solution,  then  the  activated  carbon  should  have  adsorbed  900  micro- 
moles per  gram  of  OPP  at  equilibrium.  Subsequent  addition  of  DNOSBP  at  the 
30  micromolar  level  would  result  in  no  displacement  of  the  OPP  (preferentially 
adsorbed  solute'  and  no  adsorption  of  DNOSBP  since  the  adsorption  is  assumed 
to  be  irreversible. 

EXPERIMENTA L COLUWR  PT'JVIES 

A discussion  of  twenty  experimental  studies  follows.  It  has  been  divided 
into  three  sections;  (1)  fluidized-bed  studies,  (2)  packed-bed  studies,  and  (3) 
CSTR  studies.  At  the  end  of  each  section  the  results  discussed  within  that 
section  are  summarized  in  a tabular  format  which  includes  actual  solid- 
phase  bisolute  data  as  well  as  values  predicted  using  three  models  of  multi- 
solute adsorption  equilibria.  Single-solute  data  obtained  while  the  activated 
carbon  was  being  presaturated  for  each  experimental  run  is  presented  for 
reference  in  Appendix  D of  this  report. 

Fluidized-Bed  Studies  (2, 4-DCP/DNOSBP) 

. . . r'olar  Ratio  of  l.O:  0.913  DNOSBP  to  2,4-DCP 


DNOSBP  was  the  first  solute  applied  to  the  adsorber  bed  at  an  average 
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influent  concentration  of  79.56  pM.  After  saturation,  2,4-DCP  was  added  to 
the  influent  at  an  average  concentration  of  69.37  uM,  while  the  DNOSBP 
feed  was  maintained  at  the  same  concentration.  Figure  VI I -4  and  VII-5  show 
the  solid-phase  and  fluid-phase  concentrations  versus  time,  respectively, 
for  this  run.  Table  VII-1  summarizes  actual  and  predicted  solid-phase 
concentration  equilibrium  values.  Final  equilibrium  values  for  the  run  were 
800  and  400  micromoles  per  gram  for  DNOSBP  and  2,4-DCP, respectively.  It 
is  noteworthy  that  while  the  bisolute  solution  was  fed  to  the  adsorber, 
approximately  14.9  percent  of  the  adsorbed  DNOSBP  was  desorbed  due  to  the 
addition  of  the  2,4-DCP.  The  peak  solution-phase  concentration  of  DNOSBP 
noted  after  introduction  of  the  bisolute  solution  was  approximately  107  micro- 
molar ( C/CQ  = 1.34)  and  occurred  1 hour  into  the  run.  No  re-adsorption 
of  the  DNOSBP  was  observed  after  the  initial  desorption  occurred. 

The  Langmuir  models  used  for  prediction  of  final  solid-phase  equili- 
brium values  gave  unrealistic  values  as  shown  in  Table  VII-1.  The  graphical 
irreversible  procedure  predicted  values  that  are  smaller  by  20  percent  than 
the  experimental  equilibrium  values  obtained  for  DNOSBP. 

...  Molar  Ratio  of  1,0:1.53  DNOSBP  to  2.4 -DCP 

After  pre-saturation  of  the  fl ui di zed-bed  with  79.17  uM  DNOSBP,  120.85 
urn  2,4-DCP  was  introduced  to  the  influent.  Figure  VI 1-6  shows  the  solid- 
phase  concentration  profile  for  the  run.  Figure  VII-7  indicates  the  fluid- 
phase  concentration  profile  for  the  run.  Table  VII-1  summarizes  actual 
experimental  and  predicted  solid-phase  concentration  equilibrium  values.  The 
actual  experimental  solid  phase  equilibria  values  were  520  and  900  micromoles 
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per  gram  for  DNOSBP  and  2,4-DCP,  respectively. 

After  application  of  the  2,4-DCP  it  can  be  observed  that  desorption  of 


FIGURE  VI 1-4:  SOLID  PHASE  CONCENTRATION  PROFILES  FOR  A FLUIDIZED-BED  SYSTEM 

TO  WHICH  DNOSBP  AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0:0.93  AFTER 
THE  ACTIVATED  CARBON  HAD  BEEN  PRESATURATED  WITH  DNOSBP  (79.56  uM). 


FIGURE  VI 1-5 : FLUID-PHASE  CONCENTRATION  PROFILES  FOR  A FLU  ID  I ZED-BED  SYSTEM  TO  WHICH 

DNOSBP  AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0:0.931  AFTER  THE  ACTIVATED  CARBON 
HAD  BEEN  PRESATURATED  WITH  DNOSBP  (79.56  yM). 


FIGURE  VII-6:  SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A FLU  I DI ZED-BED  SYSTEM  TO  WHICH 

DNOSBP  AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0:1.53  AFTER  THE  ACTIVATED  CARBON 

HAD  BEEN  PRESATURATED  WITH  DNOSBP  (80.0  yM) . 


FIGURE  VI 1-7:  FLUID-PHASE  CONCENTRATION  PROFILES  FOR  A FLU  I DI ZED-BED  SYSTEM  TO  WHICH 

DNOSBP  AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0:1.53  AFTER  THE  ACTIVhTED  CARBON 
HAD  BEEN  PRESATURATED  WITH  DNOSBP  (80.0  pM). 


DNOSBP  began  as  adsorption  of  2,4-DCP  occurred.  Desorption  of  DNOSBP  was  equal  to 
approximately  19.2  percent  of  the  single  solute  saturation  value.  The  peak 
solution-phase  concentration  of  DNOSBP  which  occurred  was  96  micromolar 
(C/C  = 1.2)  at  a bisolute  molar  ratio  of  1.0:1.53  DNOSBP  to  2,4-DCP. 

This  peak  occurred  at  approximately  80  hours  into  the  run.  Moreover,  during 
the  run  adsorption  and  desorption  was  observed  for  DNOSBP  although  desorp- 
tion predominated  at  final  equilibrium.  At  160  hours  after  introduction  of 
the  bi solute  mixture,  DNOSBP  adsorbed  to  a level  8 percent  greater  than  the 
presaturati on  level.  The  final  solid-phase  equilibrium  value,  however,  showed 
a net  desorption  at  steady-state. 

As  in  the  previous  case  the  competitive  and  semi-competi ti ve  Langmuir 
models  predicted  values  that  were  100  percent  greater  than  the  actual  experi- 
mental solid-phase  equilibrium  values  obtained  for  DNOSBP.  The  graphical 
irreversible  multi-solute  procedure  over-estimated  the  actual  experimental 
solid-phase  equilibrium  values  by  24  percent.  The  Langmuir  and  graphical 
predictions  for  2,4-DCP  did  not  agree  with  the  actual  values. 

...  Molar  Ratio  of  1.0:1.007  2,4-DCP  to  DNOSBP 

In  this  run  the  fl uidized-bed  contactor  was  first  equilibrated  with  80  pm 
2,4-DCP.  Thereafter,  80.6  pM  DNOSBP  was  added  as  a step  input  to  the  column. 
Figures  VII-8  and  VII-9  shows  the  solid-phase  and  fluid-phase  concentration 
profiles,  respectively,  while  Table  VI 1-1  gives  a summary  of  the  actual  and 
predicted  solid-phase  equilibrium  values.  Final  actual  solid-phase  equilibrium 
values  were  25  and  2040  micromoles  per  gram  for  DNOSBP  and  2,4-DCP,  respectively. 

The  solid-phase  concentration  profile  for  this  run  shows  that  no  desorp- 
tion of  2,4-DCP  occurred  after  application  of  DNOSBP.  This  may  partly  be 
explained  by  the  single-solute  isotherms  which  show  2,4-DCP  to  adsorb  much 
more  strongly  than  DNOSBP.  The  competitive  and  semi -competitive  Langmuir 
models  in  this  case  over-estimate  the  actual  solid-phase  equilibrium  value  for 


FIGURE  VI 1-9 : FLUID-PHASE  CONCENTRATION  PROFILES  FOR  A FLUI DIZED-BED  SYSTEM  TO  WHIC 
DNOSBP  AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1 007-1  0 AFTER  THE  ACTIVATED 
CARBON  HAD  BEEN  PRESATURATED  WITH  2,4-DCP  (80. OjjM)  . 
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the  adsorption  of  DNOSBP  by  a factor  of  ten  while  the  graphical  procedure 
predicts  no  adsorption  of  DNOSBP. 

For  2,4-DCP  the  competitive  Langmuir  model  under-estimated  the  experimen- 
tal solid-phase  equilibrium  value  by  24  percent  while  the  Semi-Competitive 
Langmuir  model  provided  an  over-estimation  by  10.6  percent.  The  graphical 
technique  indicated  an  over-estimate  of  the  actual  2,4-DCP  solid-phase 
equilibrium  concentration  of  27  percent. 

. ..  Molav  Ratio  of  1.0:2.086  2,4-DCP  to  DNOSBP 

Again,  2,4-DCP  at  80  pM  was  applied  as  the  initial  solute  to  a fluidized- 
bed  contactor.  Upon  presaturation , a step  input  of  166.86  uM  DNOSBP  was  applied 
to  the  contactor.  Figure  VII-10  shows  the  actual  solid-phase  concentration 
profile  for  this  run  after  presaturation  while  Figure  VII-11  shows  the  fluid- 
phase  concentration  profile  for  the  run.  Table  VII-1  gives  actual  and  predicted 
information  for  solid-phase  equilibria.  Final  experimental  solid-phase 
equilibrium  results  were  found  to  53  and  2130  micromoles  per  gram  for  DNOSBP 
and  2,4-DCP,  respectively. 

This  run  was  similar  in  response  to  the  previous  one  in  which  no 
desorption  of  2,4-DCP  occurred  upon  application  of  DNOSBP.  Both  Langmuir 
models  over-estimated  the  actual  solid-phase  equilibrium  value  by  almost  600 
percent.  The  graphical  technique  again  predicted  no  adsorption  of  DNOSBP. 

The  competitive  Langmuir  model  under-estimated  the  actual  solid-phase 
equilibrium  concentration  for  2,4-DCP  by  49  percent,  while  the  semi-competitive 
Langmuir  model  gave  an  over-estimate  of  only  1 percent  of  the  actual  value. 

The  graphical  procedure  over-estimated  the  actual  experimental  solid-phase 
equilibrium  value  for  2,4-DCP  by  22  percent.  This  procedure  gave  the  best 
overall  prediction  for  both  compounds  for  this  run. 


Fluidized-Bed  Studies  (OPP/DNOSBP) 


Molar  Ratio  of  1.0: 2.0  OPP  to  DNOSBP 


FIGURE  VI 1-11:  FLUID-PHASE  CONCENTRATION  PROFILES  FOR  A FLUIDIZED-BED  SYSTEM  TO  WHICH 

DNOSBP  AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  2.086:1.0  AFTER  THE  ACTIVATED 
CARBON  HAD  BEEN  PRESATURATED  WITH  2,4-DCP  (80  uM). 
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After  the  bed  of  activated  carbon  was  saturated  with  DNOSBP  at  the  80 
micromolar  solution-phase  concentration  level,  the  influent  concentration  of 
OPP  was  stepped  up  from  0 to  40  micromolar  while  the  DNOSBP  concentration  was 
maintained  constant.  Figure  V 1 1 - 1 2 shows  the  solid-phase  concentration  of 
OPP  and  DNOSBP  versus  time  while  Figure  VII-13  gives  the  corresponding 
solution-phase  concentrations.  Table  VII-1  contains  a summary  of  the  experi- 
mental and  theoretical  equilibrium  values.  Steady-state  (equilibrium) 
adsorption  values  were  observed  to  be  690  and  400  micromoles  per  gram  for 
DNOSBP  and  OPP,  respectively.  It  should  be  noted  only  a small  amount  of 
DNOSBP,  less  than  2 percent,  was  desorbed  at  steady-state.  However,  as  the  solid- 
phase  concentration  versus  time  profile  for  DNOSBP  shows,  a greater  quantity 
of  DNOSBP  (about  6.5%)  had  actually  desorbed  during  the  run  but  was  subsequently 
re-adsorbed.  The  C/CQ  versus  time  plot.  Figure  VII-13,  clearly  shows  this 
phenomenon.  This  desorption/re-adsorption  phenomenon  was  also  observed  in  a 
similar  run  in  which  the  OPP  concentration  was  80  micromolar. 

The  competitive  Langmuir  model  most  closely  predicted  the  final  OPP  solid- 
phase  concentration.  The  predicted  value  was  only  19  percent  above  the  experimental 
value.  Both  the  semi -competi ti ve  Langmuir  and  the  ideal  solution  theory  over- 
estimated the  adsorption  capacity  by  51  and  55  percent,  respectively.  The  graph- 
ical solution  method  under-estimated  the  capacity  by  the  same  amount  (50  percent). 

As  expected,  the  multi -solute  models  were  conservative  in  predicting 


equilibrium  values  for  DNOSBP.  The  competitive  Langmuir  and  semi -competitive 
Langmuir  models  were  in  error  by  50  percent.  The  ideal  solution  theory  model 
gave  an  even  more  conservative  value  and  was  in  error  by  65  percent. 


...  Molar'  Ratio  of  1.0: 1.0  OPP  to  DNOSBP 

After  equilibration  with  80  micromolar  DNOSBP,  a 80  micromolar  0PP/80 
micromolar  DNOSBP  solution  was  introduced  to  the  system.  Figures  V 1 1 - 1 4 and 
VI 1-15  show  the  solid-phase  and  solution-phase  profiles  for  OPP  and  DNOSBP, 


FIGURE  VI I -12:  SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A FLU  I DI ZED-BED  SYSTEM  TO 

WHICH  DNOSBP  AND  OPP  WERE  INTRODUCED  AT  A MUlAR  RATIO  OF  2. 0:1.0  AFTER  THE  ACTIVATED 
CARBON  HAD  BEEN  PRESATURATED  WITH  DNOSBP  (80.0yM). 
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FIGURE  V 1 1 - 1 4 : SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A FLU  I D I ZED-BED  SYSTEM  TO  WHICH 
DNOSBP  AND  OPP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1. 0:1.0  AFTER  THE  ACTIVATED  CARBON 
HAD  BEEN  PRESATURATED  WITH  DNOSBP  (80  yM) . 


respectively.  Table  VII-1  summarizes  the  comparison  with  multi-solute  equilibira 
models.  The  final  values  are  660  and  488  micromoles  per  gram  for  DNOSBP  and 
OPP,  respectively. 

In  this  as  in  the  previous  run,  a desorption-adsorption  sequence  was 
noted  for  DNOSBP.  For  this  run,  the  equilibrium  solid-phase  concentration  of 
DNSOBP  indicated  that  almost  6 percent  was  desorbed,  although  the  maximum 
desorption  occurred  at  approximately  24  hours  after  the  step  input  had  been 
initiated,  and  represented  9 percent  of  the  original  (single-solute)  capacity 
of  700  micromoles  per  gram. 

As  in  the  previous  case,  the  competitive  Langmuir  model  gave  the  best 
approximation  of  the  equilibrium  solid-phase  OPP  concentration.  It  was  in  error 
by  34  percent.  The  semi -competitive  Langmuir  and  ideal  solution  theories  predic- 
ted values  greater  than  the  actual  measured  value  by  57  and  47  percent,  respec- 
tively. The  graphical  procedure  once  again  predicted  a conservative  (small) 
value  but  was  markedly  closer  to  the  actual  value  in  this  case. 

As  above,  the  equilibrium  models  for  solid-phase  DNOSBP  concentration 
were  extremely  conservati ve.  The  competitive  Langmuir  and  semi -competitive 
Langmuir  models  predicted  values  only  36  percent  of  the  experimental ly  measured 
values.  The  ideal  solution  theory  was  less  precise,  giving  32  percent  of  the 
actual  value. 

...Molar  Ratio  of  1. 8:1.0  OPP  to  DNOSBP 

After  equilibration  of  the  activated  carbon  with  an  80  micromolar  DNOSBP 
solution,  a step  input  of  145  micromolar  OPP  was  introduced  to  the  column.  Fig- 
ures VI 1-16  and  V 1 1 - 1 7 show  the  solid  and  solution-phase  concentrations  versus 
time  for  this  run.  Table  VII-1  summarized  comparison  of  the  actual  data  with 
multi-solute  equilibria  predictions.  Final  equilibrium  values  were  545  micro- 
moles per  gram  for  DNOSBP  and  500  for  OPP. 

One  observed  difference  between  this  run  and  the  two  previous  runs  was 


DNOSBP 


FIGURE  VII-16:  SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A FLU I D I ZED-B ED  SYSTEM  TO  WHICH 

DNOSBP  AND  OPP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0: 1.8  AFTER  THE  ACTIVATED  CARBON 
HAD  BEEN  PRESATURATED  WITH  DNOSBP  (80  uM) . 
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FIGURE  V 1 1 - 1 7 : FLUID-PHASE  CONCENTRATION  PROFILES  FOR  A FLUIDIZED-BED  SYSTEM  TO  WHICH 

DNOSBP  AND  OPP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0: 1.8  AFTER  THE  ACTIVATED  CARBON 
HAD  BEEN  PRESATURATED  WITH  DNOSBP  (80  yM). 
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that  no  adsorption  of  DNOSBP  occurred  after  the  initial  desorption.  Also, 
the  amount  of  DNOSBP  desorbed  was  substantially  higher  than  for  the  two 
previous  cases.  Approximately  155  micromoles  per  gram  of  DNOSBP  were  desorbed; 
that  is,  over  22  percent  of  the  solid-phase  DNOSBP  concentration  was  displaced 
by  the  OPP. 

Another  difference  noted  between  this  run  and  the  others  was  the  period 
of  time  required  to  attain  steady-state.  The  higher  concentration  of  OPP  in 
the  influent  of  this  run  provided  a larger  driving  force  which,  in  turn,  re- 
sulted in  a more  rapid  approach  to  steady-state  conditions.  This  run  required 
approximately  100  hours  to  reach  equilibrium  while  the  three  previous  runs 
required  140  to  260  hours. 

Following  the  pattern  of  the  three  runs  described  previously,  the 
competitive  Langmuir,  semi-comoeti ti ve  Langmuir,  and  ideal  solution  theory 
predicted  values  for  solid-phase  OPP  concentrations  in  excess  of  those  experi- 
mentally observed.  Also,  the  graphical  approach  gave  values  less  than 
those  obtained  experimentally . 

Predicted  equilibrium  values  for  DNOSBP  again  were  not  accurate. 

Competitive  Langmuir  and  semi -competitive  Langmuir  predictions  were  in  error  by 
70  percent  while  the  ideal  solution  theory  prediction  was  in  error  by  77  percent. 
...Molar  Ratios  of  1.  0:2.0,  1. 0:1.0,  and  1.  8:1.0  DNOSBP  to  OPP 

Three  runs  in  which  OPP  was  applied  as  a single-solute  at  a '-oncentration 
of  80  micromolar  were  conducted.  After  presaturation,  DNOSBP  was  applied  at 
concentrations  of  40,  80,  and  145  micromolar  in  addition  to  maintaining  the  OPP 
influent  concentration  at  80  micromolar.  The  solid  and  solution-phase  concen- 
tration versus  time  profiles  for  the  1. 0:1.0  molar  ratio  case  are  shown  in  Fig- 
ure V 1 1 - 1 8 and  V 1 1 - 1 9 . Profiles  for  the  other  molar  ratio  cases  are  essentially 
identical.  Table  VII-1  summarizes  the  experimental  results  and  the  values 
predicted  by  the  multi-solute  equilibrium  models.  The  graphical  solution 
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FIGURE  VI 1-18:  SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A FLU  I D I ZED-BED  SYSTEM  TO  WHICH 
DNOSBP  AND  OPP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0: 1.0  AFTER  THE  ACTIVATED  CARBON 
HAD  BEEN  PRESATURATED  WITH  OPP  (80  uM). 


FIGURE  VI 1-19:  FLUID-PHASE  CONCENTRATION  PROFILES  FOR  A FLU  I D I ZED-BED  SYSTEM  TO  WHICH  DNOSBP 

AND  OPP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0: 1.0  AFTER  THE  ACTIVATED  CARBON  HAD  BEEN 
PRESATURATED  WITH  OPP  (80  uM). 
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is  applicable  but  predicts  no  uptake  of  DNOSBP  as  previously  discussed. 

The  solid-phase  concentrations  remained  almost  unchanged  in  the  three 
runs.  No  perceptible  quantity  of  OPP  was  desorbed  nor  was  any  additional 
quantity  of  DNOSBP  adsorbed.  The  equilibrium  solid-phase  values  for 
DNOSBP  were  6,8,  and  15  micromoles  per  gram  for  the  40,  80,  and  145  micro- 
molar influent  concentrations  of  DNOSBP,  respectively.  Equilibrium  was 
attained  in  less  than  ten  hours. 

The  competitive  Langmuir,  semi-competit‘ive  Langmuir,  and  the  ideal 
solution  theories  all  gave  conservative  estimates  for  the  solid-phase  con- 
centration of  OPP.  It  should  be  noted  that  for  these  three  runs  the 
equilibrium  solid-phase  concentration  was  the  same;  1050  micromoles  per  gram. 

As  expected,  the  predicted  values  were  more  accurate  for  the  lower  concentra- 
tions of  DNOSBP.  Errors  ranged  from  15  to  25  percent.  At  the  145  micromolar 
concentration  of  DNOSBP,  the  error  ranged  from  36  to  50  percent.  Predictions 
for  the  DNOSBP  solid-phase  concentration  were  also  imprecise.  The  multi- 
solute models  predicted  values  much  in  excess  of  the  experimental  quantities. 
Theoretical  estimates  ranged  form  13  to  23  times  the  amount  observed. 

The  graphical  irreversible  model,  however,  gave  good  estimates  of  the 
steady-state,  solid-phase  adsorption  values  that  were  experimentally  measured. 
This  indicates  that  the  assumption  of  irreversibility  is  applicable  to  the 
preferentially  adsorbed  solute. 

. . . Molar  Ratio  of  O.O: 1.0  DNOSBP  to  OPP 

This  run  was  unique  in  that  both  compounds  were  applied  as  single  solutes. 
First,  the  carbon  in  the  column  was  equilibrated  withan  80  micromolar  solution 
of  DNOSBP.  After  equilibration,  feed  of  the  DNOSBP  solution  was  replaced 
with  a feed  solution  which  had  a concentration  of  OPP  of  80  micromolar. 

Figure  VII -20  shows  the  solid-phase  concentration  of  OPP  and  DNOSBP  plotted 
versus  time.  Table  VII-1  summarizes  the  predicted  and  the  experimental  values. 
Final  equilibrium  values  were:  500  and  570  micromoles  per  gram  for  DNOSBP 


FIGURE  VI 1-20 : SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A FLUIDIZED-BED  SYSTEM 
TO  WHICH  DNOSBP  AND  OPP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  0.0: 1.0  AFTER  THE 
ACTIVATED  CARBON  HAD  BEEN  PRESATURATED  WITH  DNOSBP  (80  pM). 
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and  OPP,  respectively. 

Approximately  200  micromoles  per  gram  of  DNOSBP  were  desorbed  during  this 
run  (30  percent  of  the  DNOSBP  originally  adsorbed).  This  illustrates  the 
inapplicability  of  the  multi-solute  models  for  which  complete  reversibility 
is  assumed.  The  competitive  Langmuir  and  semi -competitive  Langmuir  models, 
for  example,  predict  that  all  of  the  DNOSBP  would  be  desorbed.  Further, 
they  predict  a higher  solid-phase  concentration  of  OPP,  980  micromoles  per 
gram,  (approximately  70  percent  more  than  was  experimentally  observed). 

The  graphical  solution  is  consistent  in  that  it  gives  conservative 
estimates  of  the  amount  of  OPP  adsorbed.  As  expected,  the  irreversible 
graphical  method  gives  an  unrealistic  low  value  which  was  equal  to  48  percent 
of  the  actual  value. 
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Packed-Bed  Studies  ( 4-DCP/DNOSBP) 

...  Molar  Ratio  of  1.0:0.92  DNOSBP  to  2,4-DCP 

A packed-bed  column  of  the  granular  activated  carbon  was  initially  saturated 
with  81.31  pM  DNOSBP  after  which  the  solute  2,4-DCP  was  introduced  to  the 
influent  at  a concentration  of  78.46  micromolar.  Figures  VI 1-21  and  V 1 1-22 
show  the  solid-phase  and  fluid-phase  time  dependent  profiles,  respectively. 

Table  VI I -2  summarizes  the  actual  experimental  solid-phase  equilibrium  values 
and  the  predicted  solid-phase  concentrations.  Experimental  solid-phase 
equilibrium  concentrations  were  found  to  be  810  and  355  micromoles  per 
gram  for  DNOSBP/and  2,4-DCP,  respectively. 

This  run,  as  for  the  fluidized  contactors,  exhibited  desorption  of  DNOSBP. 
Application  of  2,4-DCP  as  the  second  solute  caused  a 9 percent  desorption  of 
the  less  strongly  adsorbed  DNOSBP.  The  maximum  solution-phase  concentration 
of  DNOSBP  observed  due  to  desorption  was  105.4  micromolar  ( C/C Q = 1.33)  at 
a bisolute  molar  ratio  of  1.0:0.92  DNOSBP  to  2,4-DCP.  This  peak  occurred 
at  approximately  1 hour  into  the  bisolute  segment  of  the  run. 
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E VI 1-21 : SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A PACKED-BED  SYSTEM  TO  WHICH  DNOSBP 

,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0:0.92  AFTER  THE  ACTIVATED  CARBON  HAD 
PRESATURATED  WITH  DNOSBP  (81.31  yM). 
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The  competitive  and  semi -competitive  Langmuir  models  both  under- 
estimate the  solid-phase  equilibrium  value  of  DNOSBP  by  210  percent  and  over- 
estimated the  2,4-DCP  value  by  an  even  greater  error.  The  graphical  method 
predicted  values  of  80  and  535  percent  of  the  experimental  solid-phase 
equilibrium  values  obtained  for  DNOSBP  and  2,4-DCP,  respectively. 

. . . Molar  Ratio  of  1.0:1.56  DNOSBP  to  2,4-DCP 

This  run,  conducted  in  a packed-bed  contactor,  involved  equilibrating  the 
system  with  78.52  yM  DNOSBP  and  then  adding  to  the  feed  solution,  as  a step 
input,  122.76  uM  2,4-DCP.  Figure  VII -23  shows  the  solid-phase  concentration 
profiles  for  the  bisolute  portion  of  the  run.  Figure  VII -24  gives  the  solution- 
phase  concentration  profile  for  the  run  while  Table  VII-2  summarizes  actual 
experimental  and  predicted  solid-phase  equilibrium  values.  Experimental 
solid-phase  equilibrium  concentrations  were  determined  to  be  630  and  960 
micromoles  per  gram  for  DNOSBP  and  2,4-DCP,  respectively. 

As  in  the  previous  cases  of  this  sequence  of  solute  application,  a 


desorption-sorption  phenomenon  was  observed  for  DNOSBP.  At  steady-state  a 
net  desorption  of  DNOSBP  which  was  24%  less  than  the  single-solute  solid- 


phase  saturation  value  was  observed.  The  peak  fluid-phase  effluent  concentra- 
tion of  DNOSBP  due  to  desorption  was  98  micromolar  (C/C0  = 1.24)  at  a bisolute 
molar  ratio  of  1.0:1.56  DNOSBP  to  2,4-DCP.  This  peak  occurred  at  approxi- 
mately 6 hours  into  the  bisolute  portion  of  the  run. 

For  this  run  the  competitive  and  semi -competi ti ve  Langmuir  models 


predicted  values  that  were  only  30%  of  the  measured  solid-phase  equilibrium 
value  obtained  for  DNOSBP.  The  graphical  procedure  predicted  a value  which 
was  only  2 percent  higher  than  the  actual  solid-phase  concentration  obtained 


for  DNOSBP. 


Both  of  the  Langmuir  models  and  the  graphical  technique  gave  predictions 
that  were  in  error  in  excess  of  100  percent  of  the  measured  value  for  2,4-DCP. 


FIGURE  VI 1-23:  SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A PACKED-BED  SYSTEM  TO  WHICH 
DNOSBP  AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0:1.56  AFTER  THE  ACTIVATED 
CARBON  HAD  BEEN  PRESATURATED  WITH  DNOSBP  {78.52  yM). 
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FIGURE  VI 1-24:  FLUID-PHASE  CONCENTRATION  PROFILES  FOR  A PACKED-BED  SYSTEM  TO  WHICH 
DNOSBP  AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0:1.56  AFTER  THE  ACTIVATED 
CARBON  HAD  BEEN  PRESATURATED  WITH  DNOSBP  (78.52  yM). 
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...Molar  Ratio  of  1.0: 1.0  2,4-DCP  to  DNOSBP 

A packed-bed  contactor  was  brought  to  equilibrium  with  an  influent 
containing  2,4-DCP  at  a concentration  of  80  uM.  A step  input  of  80.75  pM 
DNOSBP  was  added  to  the  influent  solution  thereafter.  Figures  VII-2  and 
V 1 1 -26  show  the  solid-phase  and  fluid-phase  concentration  profiles  for  the 
bisolute  portion  of  the  run.  Table  VII-2  summarizes  experimental  and  predicted 
solid-phase  equilibrium  values.  Measured  experimental  solid-phase  equilibrium 
concentrations  were  found  to  be  29  and  1950  micromoles  per  gram  for  DNOSBP 
and  2,4-DCP,  respectively. 

Prediction  of  the  solid-phase  equilibrium  concentration  for  DNOSBP 
was  over-estimated  by  the  competitive  and  semi-competiti ve  Langmuir  models. 

The  graphical  method  predicted  no  adsorption  of  DNOSBP. 

The  model  predictions  for  2,4-DCP  were  in  much  closer  agreement  with 
the  measured  solid-phase  equilibrium  values.  The  competitive  Langmuir 
model  predicted  a value  that  was  79  percent  of  the  actual  value  while  the 
semi -competi ti ve  Langmuir  model  predicted  a saturation  value  that  was  16 
percent  higher  than  experimentally  measured.  The  graphical  procedure  gave  a 
value  that  was  33  percent  higher  than  observed. 

...Molar  Ratio  of  1.0: 2.0  2,4-DCP  to  DNOSBP 

After  pre-saturation  of  the  packed-bed  contactor  with  a 80  pM  solution 
of  2,4-DCP,  a step  input  of  160.98  pM  DNOSBP  was  added  to  the  influent.  Figure 
VI 1-27  shows  the  solid-phase  concentration  profiles  for  the  run  while  Figure 
VII -28  illustrates  the  fluid-phase  concentration  profiles.  Table  VII-2 
summarizes  the  actual  and  predicted  solid-phase  equilibrium  concentrations. 
Actual  steady-state  solid-phase  concentrations  for  the  run  were  calculated  to 
be  90  and  2110  micromoles  per  gram  for  DNOSBP  and  2.4-DCP,  respectively. 

Sliqht  disportion  of  2,4-DCP  occurred  upon  addition  of  the  DNOSBP, 
but  quickly  recovered  to  its  origianl  (single-solute)  equilibrium  value. 

The  two  predictive  models  again  failed  to  accurately  predict  the  final 
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FIGURE  V 1 1 - 2 5 : SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A PACKED-BED  SYSTEM  TO  WHICH 
DNOSBP  AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1. 0:1.0  AFTER  THE  ACTIVATED 
CARBON  HAD  BEEN  PRESATURATED  WITH  2,4-DCP  (80  yM). 


FIGURE  VII- 26 : FLUID-PHASE  CONCENTRATION  PROFILES  FOR  A PACKED-BED  SYSTEM  TO  WHICH 
DNOSBP  AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1. 0:1.0  AFTER  THE  ACTIVATED 
CARBON  HAD  BEEN  PRESATURATED  WITH  2,4-DCP  (80  uM). 
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FIGURE  VI 1-27:  SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A PACKED-BED  SYSTEM  TO  WHICH 
DNOSBP  AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  2. 0:1.0  AFTER  THE  ACTIVATED 
CARBON  HAD  BEEN  PRESATURATED  WITH  2,4-DCP  (80  uM). 


FIGURE  VI 1-28:  FLUID-PHASE  CONCENTRATION  PROFILES 
DNOSBP  AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATH 
CARBON  HAD  BEEN  PRESATURATED  WITH  2,4-DCP  (80  yM) . 
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solid-phase  concnetration  of  DNOSBP.  The  graphic  approach  predicted  no  adsorp- 
tion of  DNOSBP. 

The  models,  however,  gave  better  results  for  2,4-DCP.  The  competitive 
Langmuir  model  predicted  a value  that  was  53  percent  of  the  actual  solid- 
phase  equilibrium  value  while  the  semi -competiti ve  model  over-estimated  the 
steady-state  adsorption  value  by  19  percent.  The  graphic  technique  provided 
an  estimate  that  was  too  high  (23  percent). 
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CSTR  Studies  (2,4-DCP/DNOSBP) 

...Molar  Ratio  of  1.0:0.96  DNOSBP  to  2,4-DCP 

A CSTR  contactor  was  brought  to  equilibrium  with  a 82  pM  solution  of 
DNOSBP.  Subsequently,  a step  input  of  79  pM  2,4-DCP  was  added  to  the  influent. 
Figures  VI 1-29  and  V 1 1 - 30  show  the  solid  and  fluid-phase  concentration  pro- 
files for  the  bisolute  portion  of  this  run.  Table  V 1 1 - 3 summarizes  the  experi- 
mental and  predicted  solid-phase  equilibrium  results.  Actual  solid-phase 
equilibrium  values  were  determined  to  be  910  and  550  micromoles  per  gram  for 
DNOSBP  and  2,4-DCP,  respectively. 

As  observed  in  previous  runs  desorption  of  DNOSBP  occurred  upon  the  addi- 
tion of  2,4-DCP  to  the  influent  stream.  Steady-state  values  showed  that  11 
percent  of  the  original  single-solute  solid-phase  saturation  concentration 
of  DNOSBP  was  desorbed.  The  peak  solution-phase  concentration  of  DNOSBP  was 
noted  to  be  approximately  112  micromolar  (C/CQ  = 1.45).  This  peak  occurred 
at  150  hours  into  the  bisolute  portion  of  the  run. 

The  competitive  and  semi -competitive  Langmuir  models  did  not  provide 
realistic  predictions  as  both  predicted  concentrations  for  DNOSBP  that  were 
only  28  percent  of  the  measured  solid-phase  equilibrium  concentration.  The 
graphical  technique  predicted  a value  of  71  percent  of  the  actual  measured 
value  and  thus  provided  the  best  prediction. 

The  actual  adsorption  of  2,4-DCP  was  550  pM/g.  Both  the  Langmuir 
graphical  predictions  for  2,4-DCP  were  several  hundred  percent  above  this 
value. 

...Molar  Ratio  of  1,0:2,02  DNOSBP /2,4-DCP 

Upon  saturating  the  activated  carbon  in  the  CSTR  with  an  80  pM  solution  of 


DNOSBP,  a step  input  of  161.77  pM  2,4-DCP  was  applied  in  addition.  Figure  VI 1-31 
shows  the  solid-phase  concentration  profile  for  the  bisolute  portion  of  this 
run.  Figure  VI 1-32  gives  the  corresponding  fluid-phase  concentration  profile. 
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FIGURE  VI 1-31 : SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A CSTR  SYSTEM  TO  WHICH  DNOSBP 
AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0:2.02  AFTER  THE  ACTIVATED  CARBON 
HAD  BEEN  PRESATURATED  WITH  DNOSBP  (80  yM). 


FIGURE  VI 1-32:  FLUID-PHASE  CONCENTRATION  PROFILES  FOR  A CSTR  SYSTEM  TO  WHICH  DNOSBP 
AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  1.0:2.02  AFTER  THE  ACTIVATED  CARBON 
HAD  BEEN  PRESATURATED  WITH  DNOSBP  (80  yM) . 
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Table  V 1 1 - 3 summarizes  the  solid-phase  equilibrium  results  for  the  experimental 
and  predicted  results.  Measured  experimental  solid-phase  equilibrium  results 
were  found  to  be  820  and  1061  micromoles  per  gram  for  DNOSBP  and  2,4-DCP, 
respecti vely. 

As  in  the  previous  run  desorption  of  DNOSBP  was  observed.  The  final  solid- 
phase  equilibrium  value  was  27  percent  less  than  that  observed  for  single- 
solute saturation.  The  maximum  solution-phase  concentration  of  DNOSBP  attri- 
butable to  desorption  was  144  micromolar  (C/CQ  = 1.8).  This  peak  occurred 
at  approximately  80  hours  into  the  run. 

The  graphical  procedure  again  provided  the  best  prediction  for  the 
observed  DNOSBP  solid-phase  equilibrium  value.  The  Langmuir  models  gave 
predictions  that  were  only  20  percent  of  the  actual  value. 

Predictions  of  2,4-DCP  adsorption  levels  were  poor  for  all  models.  The 
competitive  Langmuir  model  gave  the  best  results;  90  percent  above  the  actual 
solid-phase  concentration.  The  graphical  procedure  predicted  a value  that 
was  too  high  by  109  percent. 

...Molar  Ratio  of  1.0:0.96  2,4-DCP  to  DNOSBP 

The  CSTR  was  equilibrated  with  an  82.21  uM  solution  of  2,4-DCP  before  DNOSBP 
was  added  at  the  70.93  yM  concentration  level.  Figures  VI 1-33  and  VII- 34 
show  the  solid  and  solution-phase  profiles,  respectively,  for  the  run. 

Table  V 1 1 - 3 summarizes  the  experimental  and  predicted  results.  Experimental 
solid-phase  concentration  values  obtained  for  this  run  were  79  and  2700  micro- 
moles per  gram  for  DNOSBP  and  2,4-DCP,  respectively. 

A small  quantity  of  DNOSBP  was  adsorbed  during  the  bisolute  portion  of  the 
run.  However,  the  2,4-DCP  beqan  to  adsorb  for  the  200  hours  after  the  step  input  of 
DNOSBP  was  imposed.  The  ultimate  solid-phase  concentration  of  2,4-DCP  attained 
the  equilibrium  listed  above.  This  phenomenon  was  not  observed  in  any  other  studies 


FIGURE  VII- 33 : SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A CSTR  SYSTEM  TO  WHICH  DNOSBP 
AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  0.96:1.0  AFTER  THE  ACTIVATED  CARBON 
HAD  BEEN  PRESATURATED  WITH  2,4-DCP  (82.21  yM). 
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and  cannot,  at  this  time,  be  explained. 

Again  the  Langmuir  models  over-estimated  the  adsorption  of  DNOSBP  by  a 
significant  amount  (248  percent  above  the  actual  solid-phase  equilibrium 
value).  The  graphical  method  predicted  no  DNOSBP  adsorption. 

The  competitive  Langmuir  model  under-estimated  the  2,4-DCP  solid-phase 
equilibrium  (58  percent  of  the  experimentally  measured  value.  The  semi- 
competitive  Langmuir  gave  a better  prediction  (84  percent  of  the  actual 
value).  The  graphic  method  gave  a good  approximation  at  96  percent  of  the 
experimental ly  determined  equilibrium  concentration. 

. . . Molar  Ratio  of  1.0:2.02  2,4-DCF  to  DKOSBP 

This  run  consisted  of  saturating  the  carbon  in  a CSTR  with  a 76.64 
pM  solution  of  2,4-DCP  and  then  adding  a step  of  input  of  DNOSBP  at  156  yM. 
Figures  VI 1-35  and  VII -36  show  the  solid  and  fluid-phase  concentration  profiles, 
respectively.  Table  VII-3  summarizes  the  experimental  and  predicted  solid- 
phase  steady-state  results.  Actual  experimentally  measured  values  for  DNOSBP 
and  2,4-DCP  were  225  and  2350  micromoles  per  gram,  respectively. 

During  this  run  2,4-DCP  was  observed  to  desorb.  This  was  not  detected 
in  any  of  the  previous  studies  where  2,4-DCP  was  the  initial  solute  applied 
to  the  carbon.  The  extent  of  desorption  approximated  3 percent  of  the  original 
solid-phase  saturation  value.  The  peak  solution-phase  concentration  of  2,4-DCP 
due  to  desorption  was  106  micromolar  ( C/C 0 = 1.4).  This  peak  occurred  at 

approximately  90  hours  into  the  run. 

The  Langmuir  models  both  over-estimated  the  adsorption  of  DNOSBP  by  62 
percent  of  the  measured  solid-phase  steady-state  concentration  value.  The 
graphical  irreversible  multi-solute  technique  predicted  no  adsorption  of 
DNOSBP. 

The  2,4-DCP  adsorption  predicted  by  the  competitive  Langmuir  model  amounted 


FIGURE  VI 1-35:  SOLID-PHASE  CONCENTRATION  PROFILES  FOR  A CSTR  SYSTEM  TO  WHICH  DNOSBP 
AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  2.02:1.0  AFTER  THE  ACTIVATED  CARBON 
HAD  BEEN  PRESATURATED  WITH  2,4-DCP  (76.64  uM). 


FIGURE  VI 1-36:  FLUID-PHASE  CONCENTRATION  PROFILES  FOR  A CSTR  SYSTEM  TO  WHICH  DNOSBP 
AND  2,4-DCP  WERE  INTRODUCED  AT  A MOLAR  RATIO  OF  2.02:1.0  AFTER  THE  ACTIVATED  CARBON 
HAD  BEEN  PRESATURATED  WITH  2,4-DCP  (76.64  pM). 
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to  only  47  percent  of  the  experimentally  measured  solid-phase  equilibrium 
value  while  the  semi-competi ti ve  model  was  in  much  closer  agreement  (91  percent 
of  the  measured  value). 

SUMMARY 

Observation  of  the  experimental  results  and  examination  of  the  solid- 
phase  and  fluid-phase  concentration  profiles  indicate  a positive  preferential 
adsorption  for  2,4-DCP  and  OPP  over  DNOSBP.  This  was  to  be  expected  as  the 
isotherms  for  the  compounds  indicate  2,4-DCP  and  OPP  adsorb  in  much  greater 
amounts  than  does  DNOSBP  when  being  applied  to  carbon  in  equal  concentrations. 
This  preferential  adsorption  can,  in  part,  explain  the  ability  of  2,4-DfF  or 
OPP  to  displace  DNOSBP  from  the  absorbent  during  a run.  This  also  explains 
the  observation  that  2,4-DCP  and  OPP,  tended  to  be  irreversibly  adsorbed  and 
did  not  desorb  as  a result  of  adding  DNOSBP  to  the  experimental  system. 

Calculations  using  the  Langmuir  multi-solute  models  and  the  graphical 
procedure,  showed  that  neither  can  be  accepted  as  logical  tools  for  predicting 
multi-solute  equilibria  under  dynamic  condition. 

I! 


[ 


CHAPTER  VIII 

MATHEMATICAL  SIMULATIONS 
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A series  of  model  simulations  were  conducted  for  the  purpose  of 
illustrating  the  design  and  operational  alternatives  one  has  for  minimizing 
the  occurrence  of  chromatographic  displacement  of  adsorbed  contaminants  from 
an  adsorber  column.  This  involved  use  of  the  mathematical  model  developed  in 
Chapter  VI  which  consists  of  a set  of  dispersed  flow  mass  balance  equations 
and  the  appropriate  kinetic  and  equilibrium  relationships. 

SIMULATION  CONDITIONS 

For  the  current  simulations  film  diffusion  was  considered  to  limit  the 
rate  of  transport  of  solute  molecules  from  the  bulk  solution  to  the  adsorbent 
surface.  Because  neither  the  Langmuir  Competitive,  Langmuir  Semi-Competi ti ve. 
Ideal  Solution  Theory,  or  Graphical  adsorption  isotherm  models  adequately 
described  the  multi -solute  competitive  adsorption  equilibria  as  shown  in 
Chapter  VII,  the  Langmuir  Competitive  Adsorption  isotherm  model  which  assumes 
complete  reversibility  was  employed  for  the  simulations  conducted  as  part  of 
this  study.  This  selection  was  made  because  the  Langmuir  Competitive  multi- 
solute adsorption  model  allows  one  to  simulate  the  most  critical  elution 
conditions  that  could  potentially  occur  in  prototype  adsorbers.  This  in 
turn,  permits  the  specification  of  a conservative  (fail  safe)  set  of  design 
and  operational  criteria. 

To  simulate  the  most  critical  elution  conditions  from  an  engineering  design 

and  operations  criteria  viewpoint,  all  simulations  were  conducted  for  the  case 

of  two  competing  solutes,  parani trophenol  (PNP)  and  parabromophenol  (PBP)  which 

were  introduced  to  the  adsorber  column  at  concentrations  of  l.OmM.  For 

simulations  the  columnar  adsorber  was  assumed  to  contain  a polymeric  adsorbent, 

2 

Duolite  A-7,  and  to  be  operated  at  a flow  rate  of  5 gpm/ft  . The  Langmuir 
equilibrium  constants  employed  for  the  simulations  are: 
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Parani trophenol  (PNP) 

Q = 1.515  millimoles  per  gram 
b = 7457  liters  per  mole 

Parabromophenol  (PBP) 

Q = 3.495  millimoles  per  gram 
b = 650  liters  per  mole 

It  is  noteworthy  that  parabromophenol  is  adsorbed  to  a greater  extent  than 
is  parani trophenol  at  ultimate  monolayer  saturation,  but  has  a lower 
energy  of  adsorption  as  evidenced  by  the  relative  b values. 

SIMULATIONS 

Six  different  simulations  were  conducted  as  part  of  this  study. 

The  first  two  were  conducted  for  the  operational  condition  in  which  the 
adsorbent  was  first  saturated  with  PBP  as  a single-solute.  After  saturation, 
both  PNP  and  PBP  were  introduced  to  the  column  at  concentrations  of  1.0  mM. 
The  two  simulations  were  conducted  for  the  packed-  and  fluidized-bed 
cases.  For  the  second  set  of  simulations  the  order  of  addition  of  the  com- 
peting compounds  was  reversed.  That  is,  PNP  was  introduced  first  as  a 
single-solute  and  the  bisolute  mixture  was  added  thereafter.  These,  again, 
were  run  for  the  packed-  and  fluidized-bed  cases.  The  final  set  of  simula- 
tions were  conducted  for  the  case  in  which  both  PNP  and  PBP  were  introduced 
to  the  column  simultaneously.  As  before,  these  simulations  were  conducted 
for  packed-  and  fluidized-bed  modes  of  operation. 

PBP  PNP 

A set  of  simulations  was  conducted,  as  noted  above,  for  the  case  in 
which  the  adsorbent  in  the  columnar  adsorber  was  presaturated  with  PBP  at 
the  l.OmM  concentration  level.  Subsequently,  a solution  which  contained 
both  PNP  and  PBP  at  the  1.0  mM  level  was  introduced  to  the  column.  The 
simulation  traces  for  both  solution  and  solid-phases  are  given  in  Figures 
VI 1 1 - 1 and  VIII-2  for  the  packed-  and  fluidized-bed  cases,  respectively. 
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It  can  be  observed  that  as  PNP  was  introduced  to  both  packed-  and 
fluidized-bed  adsorbers  a considerable  quantity  of  adsorbed  PBP  was  dis- 
placed from  the  adsorbent  and  was  eluted  into  the  effluent.  This  is  evi- 
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denced  by  the  concentration  peak  which  can  be  observed  for  PBP  in  the  solution 
phase  traces  and  by  the  decreasing  solid-phase  concentration  of  PBP.  Further 
analysis  of  these  simulations  shows  that  PNP  is  the  preferentially  adsorbed 
solute  and  that  PBP  is  competitively  adsorbed.  Comparison  of  the  packed-bed 
and  fluidized-bed  profiles  showed  only  minor  differences:  the  fluidized-bed 

solution-phase  elution  peak  had  a slightly  greater  amplitude  but  did  not  pre- 
vail as  long  as  the  elution  peak  for  the  packed-bed  case. 

PNP  — ' -PBP 

For  this  pair  of  simulations  for  packed-  and  fluidized-bed  adsorbers, 
the  adsorbed  was  first  presaturated  with  PNP  at  the  l.OmM  concentration 
level.  A solution  which  contained  both  PNP  and  PBP  at  1.0  millimole  per 
liter  concentrations  was  subsequently  introduced  to  the  column.  Figures 
VI 1 1 -3  and  VIII-4  show  the  solution-  and  solid-phase  profiles  for  this 
series  of  simulations. 

The  simulation  traces  show  that  some  PNP  is  displaced  from  the  adsorb- 
ent and  is  thereby  eluted  from  the  columnar  adsorber.  This  occurred  for 
both  the  packed-  and  fluidized-bed  cases.  Nonetheless,  it  is  extremely 
important  to  point  out  that  the  amplitude  of  the  PNP  elution  peak  was 
dampened  considerably  in  the  fluidized-bed  adsorber.  These  simulations  also 
confirm  the  prior  observation  that  PNP  is  the  preferentially  adsorbed  solute 
and  that  PBP  is  competitively  adsorbed. 

Comparison  of  Figures  VI 1 1 -3  and  VIII-4  with  Figures  VIII-1  and  VIII-2 
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show  that  the  magnitude  of  the  elution  peak  and  the  total  quantity  of 
solute  eluted  are  diminished  considerably  when  the  preferentially  adsorbed 
solute  is  applied  to  the  column  prior  to  the  competitively  adsorbed  solute. 
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FIGURE  VIII-4:  SOLUTION-  AND  SOLID-PHASE 

CONCENTRATION  PROFILES  FOR  PBP  AND  PNP  FOR 
A FLUIDIZE D-BED  ADSORBER 


119 

This  observation  has  rather  significant  implications  with  respect  to  the 
operation  of  prototype  adsorbers.  That  is,  if  possible,  production  schedules 
should  be  arranged  such  that  wastewaters  containing  the  highest  energy  ad- 
sorbing contaminants  are  introduced  first  to  a column  of  fresh  adsorbent  and 
the  wastewaters  containing  the  lowest  energy  adsorbing  contaminants  are 
introduced  last. 

PBP  & PNP 

The  final  series  of  simulations  was  conducted  for  the  case  in  which 
both  PBP  and  PNP  were  introduced  to  the  adsorber  column  simultaneously  at 
concentrations  of  l.OmM.  Figures  VII 1-5  and  VI 1 1 -6  show  the  resultant 
solution-  and  solid-phase  traces  for  the  packed-  and  fluidized-bed  cases, 
respectively. 

For  the  packed-bed  mode  of  operation,  it  is  noted  that  a small 
quantity  of  PBP  was  displaced  from  the  adsorbent  and  was  subsequently 
eluted.  This  occurred  due  to  the  relative  rates  of  adsorption  of  the  two 
contaminants  which  resulted  in  a chromatographic  separation  of  the  compounds 
within  the  column.  That  is,  PNP,  being  the  preferentially  adsorbed  solute, 
adsorbed  near  the  influent  end  of  the  column  at  the  outset  of  the  run  while 
the  PBP,  being  the  competitively  adsorbed  solute,  was  adsorbed  primarily 
near  the  effluent  end  of  the  adsorber.  As  time  proceeded  the  PNP  moved  as 
a wave  front  through  the  column  displacing  the  PBP  ahead  of  it. 

Comparing  the  results  for  the  fluidized-bed  mode  of  operation  with 
those  obtained  for  the  packed-bed  mode,  one  observes  that  essentially  no 
PNP  is  displaced  from  the  adsorbent.  This  can  be  attributed  to  the  fact  that 
the  adsorbent  solids  are  virtually  completely  mixed  in  a fluidized  columnar 
adsorbent.  This  prevents  the  chromatographic  separation  of  the  solutes 
within  the  column  and,  therefore,  eliminates  the  possibility  for  displace- 
ment and  elution. 
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FIGURE  VI 1 1-5:  SOLUTION-  AND  SOLID-PHASE 

CONCENTRATION  PROFILES  FOR  PBP  AND  PNP  FOR 
A PACKED-BED  ADSORBER 
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FIGURE  VIII -6 : SOLUTION-  AND  SOLID-PHASE 

CONCENTRATION  PROFILES  FOR  PBP  AND  PNP  FOR 
A FLUIDIZED-BED  ADSORBER 
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These  observations  again  underscore  the  conclusion  that  the  fluidized- 
bed  mode  of  operation  has  distinct  operational  advantages  over  the  packed- 
bed  method  of  sol ids-1 iquid  contact.  Moreover,  the  results  of  this  series 
of  simulations  clearly  establish  that  the  effects  of  chromatographic  dis- 
placement and  elution  can  be  minimized  if  the  competing  solutes  are  intro- 
duced to  the  adsorber  simultaneously.  This  could  be  accomplished  by  arrang- 
ing industrial  production  schedules  such  that  all  contaminants  are  discharged 
concurrently.  If  this  is  not  possible,  however,  then  provision  should  be 
made  f c * concentration  equilization  ahead  of  the  adsorption  system. 

SYNOPSIS  OF  SIMULATIONS 

Results  of  the  six  simulations  conducted  have  been  summarized  in 
Table  VI 1 1 - 1 . The  following  observations  can  be  tabulated: 

(1)  If  the  competitively  adsorbed  solute  (PBP)  is  introduced 
to  the  adsorber  first,  then  considerably  more  solute  is 
eluted  when  the  preferentially  adsorbed  solute  is  introduced 
to  the  column  than  if  the  order  of  addition  were  reversed. 

(2)  Due  to  mixing  of  the  adsorbent  solids,  the  fluidized-bed 
mode  of  operation  is  substantially  better  in  minimizing  the 
effects  of  contaminant  displacement  and  elution. 

(3)  When  all  competing  contaminants  are  introduced  to  the 
adsorber  simultaneously,  the  effects  of  displacement  and 
elution  are  minimized. 


CHAPTER  IX 


CONCLUSIONS  AND  RECOMMENDATIONS 

Neither  the  Langmuir  Competitive,  Langmuir  Semi -Competitive,  Ideal 
Solution  Theory,  or  Graphical  adsorption  isotherm  models  adequately 
described  multisolute  competitive  adsorption  equilibria. 

To  minimize  the  elution  of  adsorbed  contaminants  from  adsorbers 
the  fluidized-bed  mode  of  operation  offers  distinct  advantages 
over  packed-bed  operation. 

If  possible,  industrial  production  schedules  should  be  arranged 
such  that  all  contaminants  are  discharged  concurrently.  This  also 
serves  to  minimize  the  effects  of  chromatographic  displacement  of 
adsorbed  contaminants. 

If  the  above  (3)  is  not  possible,  then  production  schedules  shoulu 
be  arranged  such  that  wastewaters  containing  the  highest  energy 
adsorbing  contaminants  are  introduced  first  to  a column  of  fresh 
adsorbent  while  the  wastewaters  containing  the  lowest  energy  adsorb- 
ing contaminants  are  introduced  last. 

Another  design  option  which  serves  to  minimize  displacement  and 
subsequent  elution  adsorbed  contaminants  is  to  make  provision  for 
concentration  equilization  in  the  process  configuration  ahead 
of  the  columnar  adsorbers. 
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B£T  energy  tern,  liters/mole 
cross-sectional  area  of  column,  sq.  cm. 

Langrrn  ir  energy  tern,  liters/mole 
solution-phase  concentration  of  solute,  moles/ liter 
saturation  concentre t ion  of  solute  in  solution,  moles/1 itor 
pore  d i f tus  iv  i ly,  cVy  sec 

p/3 

(Kj/U)  * (v/l>rf ) , mass  transfer  factor,  dimensionless 

adjustatile  turve-f ittiny  constant 
adjustable  curve-f i tt ing  constant 
number  of  elements 

U • D /v,  Reynolds  number,  dimensionless 

v/D^,  Schmidt  number,  dimensionless 

K,  * D /D„ , Sherwood  number,  dimensionless 

ultimate  uptake  capacity  of  adsorbent,  moles/gram 

solid-phase  concentration  of  solute,  moles/gran 

initial  solid-phase  concentration  (usually  zero),  moles/gram 

separation  factor,  dimensionless 

rate  of  adsorption,  moles/gram-hour 

radial  Jistance  within  particle 

times,  hours 

solution  volumetric  flow  ra.te,  liters/hour 
U/A,  ve.ocity  of  flow,  cm/fcr 
volume  of  the  element,  liters 
axial  distance,  cm 
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void  ratio  or  porosity,  dir.ensionl ess 

1/(R  * 15  (1  - R)/-.2) 

? 

Kinematic  viscosity  of  solution,  cm  /sec 

packed-bed  density,  grams/liter 

density  of  adsorbent  particle,  gm/liter 

internal  porosity  of  adsorbent  particles,  dimensionless 


DESCRIPTION 


Calgon  Corporation's  granular  activated  carbon  products  (Filtrasorb 
200)  are  made  from  select  grades  of  coal  resulting  in  superior  hardness 
and  long  life.  Produced  under  rigidly  controlled  conditions  by  high 
temperature  steam  activation,  they  are  high  density  carbons  with  high 
surface  areas.  Their  pore  structures  have  been  carefully  controlled  for 
the  adsorption  of  both  high  and  low  molecular  weight  impurities  from 
water.  These  coal-based  carbons  wet  readily,  do  not  float,  and  their 
high  density  allows  convenient  backwashing  under  conventional  flow  rate 
conditions.  Each  is  available  in  a different  mesh  size  to  suit  specific 
engineering  design  requirements. 

PURPOSES 

Filtrasorb  200  granular  activated  carbon  products  are  designed  for 
efficient  use  in  water  filtration  equipment  for  the  removal  of  organic 
impurities  found  in  potable  and  industrial  water  supplies.  These  impuri- 
ties include  tastes,  odors,  color,  insecticides,  detergents,  phenols, 
and  other  contaminants  from  both  natural  or  industrial  sources.  Further- 
more, Filtrasorb  200  is  an  effective  filtration  media  for  removal  of 
turbidity  and  suspended  matter.  In  addition  to  the  adsorption  and  fil- 
tration functions,  they  are  capable  of  concurrently  removing  excess 
oxidants  such  as  chlorine  and  permanganate  used  in  pretreatment  applica- 
tion. 


REACTIVATION 

Granular  activated  carbon  has  been  specifically  manufactured  for 
thermal  reactivation  and  re-use.  Numerous  reactivation  installations 
have  demonstrated  the  feasibility  and  economy  of  granular  carbon 


reactivation.  The  carbon  is  handled  as  a water  slurry  between  the 


filters  and  the  furnace  area.  No  additional  labor  or  labor  skills  are 


required  to  operate  a reactivation  facility  which  is  automated  to  run 


with  minimal  surveillance 


PHYSICAL  PROPERTIES 


FILTRASORB 

200 


Total  surface  area  (N2,  BET  method)  m2/g 
Bed  density  backwashed  and  drained,  lbs/ft 
Particle  density  wetted  in  water  g/cc  . . 

Effective  size  mm  

Uniformity  coefficient  


800  - 900 
30  approx. 
1.4  - 1.5 
0.55  - 0.65 
1 . 7 or  less 


SPECIFICATIONS 


FILTRASORB 

200 


Sieve  size  U.  S.  Std.  series 
Larger  than  No.  8- -Max  % 
Larger  than  No.  14--Max  % 
Smaller  than  No.  30--Max 
Smaller  than  No.  40--Max 
Mean  particle  dia.  mm  ...  . 

Iodine  number--Min  

Abrasion  number- -Min  .... 

Ash,  Max  % 

Moisture  as  packed.  Max  % . . 


Calgon  Corporation,  Activated  Carbon  Product  Bulletin  20-1 


The  following  program  was  used  to  compute  the  coefficients  for  the 


fifth  order  polynomial  in  describing  single-solute  isotherms: 


C PROGRAM  FOR  SOLVING  SIMULTANEOUS 
C ALGEBRAIC  EQUATIONS  BY  USE  OF 
C SUBROUTINE  SIMQ 

DIMENSION  A(100),  X(10) 

C READ  THE  NUMBER  OF  EQUATIONS 
READ  (1,100)  INDEX 
100  FORMAT  (12) 

C READ  THE  COEFFICIENTS  AND  THE  X 
C ARRAY 

DO  10  J = 1,  INDEX 
IK  = INDEX  * (INDEX  - 1)  + J 
READ  (1 , 110) (A(I) , I=J , IK . INDEX), X(J) 
110  FORMAT  (6F12.0) 

10  WRITE  (3,120)  (A(I)=J, IK, INDEX) ,X(J) 
120  FORMAT  (1HO, 11F12.0) 

CALL  SIMQ(A,X, INDEX, KS) 

WRITE(3, 130)  KS 
130  FORMAT (1H,// , 3HKS=, 14) 

DO  20  1=1,  INDEX 
20  WRITE(3, 140)1, X(I) 

140  FORMAT ( 1H , 10X, 2HX( , 12 , 2H) = , E15 . 6) 

STOP 

END 


Note:  The  isotherm  is  separated  iito  zones.  For  example,  the  simul- 

taneous equations  were: 


(X/m)j  = KiCej  + 
(X/m) 2 = KiCe2  + 

(X/m) 5 = KiCe5  + 


K2Ce  i + K3Cei 
K2Ce|  + K3Ce2 

K2Ce 5 + K3CeI 


+ KuCei  + KsCe? 
+ Ki»Ce2  + KjCe2 

+ KuCes  + KsCel 


The  following  computer  program  represents  the  technique  used  for 
determining  the  spreading  pressure  from  single-solute  equilibrium  data: 


$$$ COST  I NUOUS  SYSTEM  MODELING  PROGRAM  III  V1M0  TRANSLATOR  OUTPUTSSS 

* COMPUTATION  OF  RADKE  AND  PRAUNITZ 
INITIAL 

RENAME  TIME=  CONC 
INCON  CONC  =1 . E-12 
DYNAMIC 

FUNCTION  OPP  = (0.00,0. 00), (5,1385), (10, 1545), (15, 1655), (20,1745),... 
(25,1820) , (30,1880) , (35,1938) , (40,1988) , (50,2070) , (60,2140) , . . . 

(80,2234)  , (90,2260) , (100,2280) , (110,2295) , (120, 2305) , (130, 2310) .... 
(140,2318) , (160,2320) , (130,2325) , (200 , 2325) , (220 , 2325) 

FUNCTION  DNOS  = (0 . 0, 0 . 0) , (5 , 870) , (10, 1015) , ( 15 , 1 100) , (20 , 1 1 55) , . . . 
(25,1195) , (30,1225) , ( 35 , 1255) , (40 , 1 280) , (50,1310) , (60,1340) .... 

(80,1380) , (90,1395) , (100,1405) , (110,1415) , (120,1420) , (130,1428) .... 
(140,1430) , (160,1440) , (180, 1450) ,(200,1450) ,(220,1450) 

NOC  =NLFGEN (OPP, CONC) 

NDC  =NLFGEN( DNOS, CONC) 

NOPP  = NOC/CONC 
NDNS  = NDC, CONC 
OPPT  = I NTGRL (CONC , NOPP) 

DNOST  = INTGRL (CONC, NDNS) 

TIMER  FINTIM  = 225,  PRDEL  = 1 

METHOD  RKS 

PRINT  OPPT, DNOST 

END 

STOP 


as:: 


10000 


1000 


The  method  of  solution  of  multi-solute  equilibrium  system  per 


telephonic  communication  Dr.  C.  Radke: 


FORTRAN  COMPILER  VERSION  2.3  B.l  15  FEB  71 


SUBROUTINE  NEWTON  ( NMAX,  X,  Cl,  C2  ) 

COMMON  Al,  A2,  Bl,  B7,  RN1 , RN2 
DO  1 1=1, NMAX 

CALL  FUNC  ( I,  X,  F,  DF,  Cl,  C2  ) 

XN  = X - F/DF 

S = ( ( XN  - X )**2  ) / XN**2 
X = XN 

IF  ( S - l.E-5  ) 2,2,5 

5 IF  ( X . GT.  1.0  ) X = 0.95 

IF  ( X . LT.  0.0  ) X = 0.05 

1 CONTINUE 

2 PRINT  3 

3 FORMAT  ( //  * I MOLE  FRAC  FUNCTION  ERROR  *//) 
PRINT  4,  I,  X,  F,  X 

4 FORMAT  ( 14,  3F.12.5//  ) 

RETURN 

END 


NMAX  Maximum  iterations  allowed  for  Newton  procedure  initial  guess 

X and  final  answer  for  adsorbed  phase  composition  (use  a guess 

of  X = 0.5) 

Cl,  C2  Given  mixture  concentrations  of  solutes  1 and  2. 


Al,  A2,  Constants  in  equations  H2  (32). 

Bl,  B2, 

RN1,  RN1 

FUNC  Subroutine 

t"  Vi 

I The  i iteration  of  Newton  procedure. 

Ci  C2 

F The  function  whose  root  is  desired  f ( Z i ) =0=F  !■=—  I -f  z U — — 

Z\  l-Z 

DF  The  derivative  of  that  function. 


XN  The  next  best  root. 

S Test  of  convergence.  When  S < 1 0~ 5 , X is  the  answer  you  may 

want  another  test  other  than  percentage  error. 
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FORTRAN  COMPILER  VERSION  2.3  B.l 


SUBROUTINE  FUNC  ( I,  X,  F,  DF,  Cl,  C2  ) 

COMMON  Al,  A2 , Bl,  B2,  RN1,  RN2 
Z = X 

F = Z/(A1*C1)  - (1 . -Z) / (A2*C2)  + ( (Z/C1)**RN1)/B1  - C ( 

1 C2) **RN2)/B2 

DF  = 1 . / (A1*C1)  + I . / (A2*C2)  + ( RN1* (Z** (RN1- 1 .)))/( (Cl* 
1 + ( RN2*((1.-Z** (RN2-1 . 1 ))/ ( (C2**RN2) *B2) 

RETURN 

END 


iiifidb  rtf  ini'itf 


(l.-Z)/ 
*RN1) *B1) 
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*************************************** ********************************* 
************************************************************************ 
**************  ********************************************************** 


initial 

************************************************************************ 

FIXED  I,J,< 

STOKAGE  CAUN.CBIIU 

************************************************************************ 
***  PARAMETER  INPUTS  REQUIRED  ARE  DEFINED  AS  FOLLOWS  *♦* 

» *************************  ***  ************************************  ******* 


CARBON  MASS  OF  CARBON  IN  BED  (GRAMS!  *♦* 

A EXTERNAL  TRANSFER  AREA  FOR  CARBON  (SQ  CM/GRAM ) *** 

0p  DIAMETER  OF  CARBON  PARTICLES  (CM) 

RHDPAC  PACKED  BED  DENSITY  (GRAMS/LITER)  ♦** 

EPS  PACKEO  BED  POROSITY  (DIMENSIONLESS)  *** 

QAMAX  LANGMUIR  ULTIMATE  UPTAKE  CAPACITY  FOR  A ( MOL E S/ GR AM ) *** 

QBMAX  LANGMUIR  ULTIMATE  UPTAKE  CAPACITY  FOR  B (MOLES/GRAM)*** 

BA  LANGMUIR  ENERGY  TERM  FOR  A (LITERS/MOLE) 

BB  LANGMUIR  ENERGY  TERM  FOR  B ( L I TE RS/MOL E ) **♦ 

C A ( 1 ) SOLUTI ON-PHASE  CONCENTRATION  OF  A ( MOLES/L  I TER  ) *♦* 

CB(l)  SOLUTION-PHASE  CONCENTRATION  OF  B (MOLES/LITER)  *** 

DLA  D I FFUS I V I T Y OF  SOLUTE  A (SQ  CM/SEC) 

DL3  DIFFUSIVITY  OF  SOLUTE  B ( SQ  CM/SEC)  **♦ 

AREA  CROSS-SECTIONAL  AREA  OF  COLUMN  (SO  CM)  **♦ 

j SOLUTION  VOLUMETRIC  FLOW  RATE  (LITERS/HOUR)  *** 

GAMMA  KINEMATIC  VISCOSITY  OF  WATER  (SQ  CM/SEC)  *** 

EXP  RATIOS  EXPANOED  /PACKED  BED  VOLUME  (DIMENSIONLESS)  *** 


****************************************************************** 

PARAMETER  INPUTS  FOLLOW 

****************************************************************** 


PARAM 

CARB0N=36. 79 

PARAM 

A=  124.6 

PARAM 

DP=9.S5E-2 

PARAM 

RHOPAC=253 • 

PARAM 

EPS=D. 43 

PARAM 

QAWAX=  3.4958E -3 

PARAM 

BA=65D . 

TABLE 

CAI l )=l.E-3 

TABLE 

CB(  l >= l.E-3 

PARAM 

CAQ=l.E-3 

PARAM 

CB0=l.E-3 

PARAM 

DL  A = B.  672E-6 

OARAM 

OL  8=  8. 54E-6 

PARAM 

ARE  A = 4 . 9 78 

PARAM 

U*6. 12 

PARAM 

GAMMA=9.348E- 3 

PARAM 

EXP>I l.Of 1.5) 

_ - 


********************************************************************* 

***  CALCULATED  INPUT  PARAMETERS  ARE  DEFINED  AS  FOLLOWS 

**************  ************************************* ****************** 


*** 

***  RHD  EXPANDED  BED  DENSITY  (GRAMS/  LITER) 

**«  Vp  TOTAL  VOLUME  OF  PACKED  BED  (LITERS) 

***  v TOTAL  VOLUME  OF  PACKED  BED  (LITERS) 

***  EPSI  EXPANDED  BED  POROSITY  ( D I ME  NS  I ONLC SS ) 

***  UBAR  SUPERFICIAL  VELOCITY  OF  FLOW  THROUGH  BED  (CM/SEC) 

***  NREMOD  MODIFIED  REYNOLDS  NUMBER  (DIMENSIONLESS) 

***  N SC  A SCHMIDT  NUMBER  FOR  A (DIMENSIONLESS) 

***  NSCB  SCHMIDT  NUMBER  FDR  B (DIMENSIONLESS) 

***  JD  MASS  TRANSFER  FACTOR  (DIMENSIONLESS) 

***  KF A mass  TRANSFFR  COEFFICIENT  FOR  A (CM/HOUR) 

***  <FB  MASS  TRANSFER  COEFFICIENT  FOR  8 (CM/HOUR) 

***  DUMl  DUMMY  COEFFICIENT 

***  DUM2  DUMMY  COEFFICIENT 

***  DUMB  DUMMY  COEFFICIENT 

*** 


ft***«tftft*ft***ftft*ft*ft*********ft*******ft****«*ft******* ****************** 

***  INPUT  PARAMETER  CALCULATIONS  FOLLOW 

** *********  ********************************************************** 

RHO=RHOPAC/5XP 

VP  = CARBDM/( R HOP  AC *10.  ) 

V *C  ARBON/ ( RHO * l 0 . ) 

EPSI  = (EPS*VP«-V-VP)/V 
UBAR=(U*0. 277777) /AREA 
NREMOD=(UBAR*DP)/(GAMMA*( l.-EPSI  ) ) 

NSCA=GAMMA/DL A 
NSC3=GAMMA/DLB 
JD=  L . 34*NREMDD**(  -0 . 46R  ) 

<FA=( JD*JBAR/NSCA**(2./1. ) ) *3600. 

KFB=( JD*UBAR/NSCB**( 2./ 3. ) )*3600. 

DUMl  = U/(E  PSl *V) 

0UM2= (KFA*A*RH0*1 .E-3 )/EPSI 
DUM3=(KFB*A*RHO*l.E-3)/EPSI 


***  INITIAL  CONDITIONS  FOR  SOLID  PHASE  CONCENTRATION  *** 

« ****************************************  ******************************* 


TABLE  ICOAI 2-L 1 )= 11*0. D 
TABLE  IC0B(2-ll ) = 1 1 *0  • 0 


*** 


FUNCTION  I SDB  = ( 410. « 500. ) , < 530. , l 000. ) , ( 62 8 . , 1 500 . ) , . . . 

(710. ,2000.), (783. , 250 3 . ) , ( B5 1 . , 3000. ) , ( 012. ,3500. ) , . .. 

( 0 70. ,*000. > . ( 10  75. , 5000. ) , { ll 70. ,6000. ). ( l 3* 7., 8 003.  ) , ... 

( 152  3.,  1 00  00. ),(  1700. , 12000. ) , ( 1877. ,14000. ),<2 054. ,16000.),... 
( 22  30. ,18D00. ) , ( 24D7. ,20000. ) , ( 4000. , 3B0  31 . ) , ( 6000 ., 60668 . ) ,.  .. 
( 8000. , 83  305.  ) , ( l 0000. , 105  04  2. ), ( 12  000.,  12  85  79.  ),..  . 

( 14000. , 151216. ) * ( 16000. , 173853. ), ( 18000. ,196491 .),... 

I 20000.,  21  H2B.  ) 


« i 
*** 


LINEARIZED  ISOTHERM  DATA  FOR  COMPONENT  B (PNP) 


142 


*******♦♦♦*♦***********************•*♦*****♦***♦***♦♦♦♦♦•♦♦♦♦•♦♦*♦♦***♦♦ 


)VNA*IC 

******* ***************************************************************** 

***  DIFFERENTIAL  EQUATIONS  FOR  AOSORPTION  BED  *♦* 

************************************************************************ 

qa=intgrli icqa.qadot, u ) 

QB= I NTGRL ( I CQB, QBD3T , l l ) 

PROCEDURE  QAOOT,OBOOT  = ADSOA,QB,ISOB,QAMAX,BA,KFA,KFB,A)  <j 

DO  10  1=2,11 

IF  ( QB  ( II.LT.5.E-5)  Q8INV=2.E«^ 

IF  ( QB ( I 1 • SE • 5. E-5 ) QB I NV= l /Q8 < l 1 
CBI NV=FUNGEN( I SOB , 5 , QB I N V ) 

RSLOPE  = S LOPE  ( I SOB, 3, QBl NV ) 

QBMAX  = l/<  QBI NV-( 1/BSLOPE )*CBINV ) 

BB=8SL0PE/Q3NAX 

CSTARA  = QBMAX*BB*QA(  11/ ( QAMAX*QBMAX*BA*SB-QAMAX*QB< I I *BA*BB. . . 

-QBMAX*QA( l I *BA*BB I 

CSTARB=QAMAX*BA*QB( I ) /( QAMAX*QBMAX*8A*3B-QAMAX*QB( I ) *B A*BB . . . 

-QBHAX*QA( I I *8A*BB) 

C A f I I = ( OUMl  *C  A ( t-H  *0UM2*CSTARAI  / ( DJM1  «-DUM2  I 
C B ( I )=(0UM1*CB(  I-l)«-UUM3*CSTARB) /(DUMI+0UM3) 

QADOT ( I»=XFA*A*(CAII)-CSTARA)*l.E-3 
QBDOTI I )=KFB*A*(CB(I) -CSTAR8 I *1 . E-3 

10  CONTINUE 
ENDPROCEOURE 

*******************  ************************************************** 

***  TO  ACCOUNT  c OR  MIXING  OF  CARBON  IN  FLUIDIZED  BEDS  (EXP>1.0), 

***  AVERAGE  SOLID-PHASE  CONCENTRATIONS  ARE  CALCULATED  AND  QA*  S L 

***  OB'S  ARE  REINITIALIZED  AT  THE  AVERAGE  VALUES 

********************************************************************* 

PROCEDURE  QAT,QBT=SORBED(QA,QB) 

QAT0T=0. 

3BT3T  = 0. 

DO  LI  J=?,ll 
QAT3T=QAT0T*QA( J I 

11  QBTOT=QBTOT*QBI J) 

QAT=3AT0T/10. 

QBT  = QBT3T / 10 . 

IF  (EXP. EQ. 1.0)  GO  TO  13 
DO  IP  K » 2 » l l 
3A ( < I =3 A T 

12  QB( < 1 =QBT 

13  CONTINUE 
ENDPROCEOURE 
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ft*********************************************************************** 

TERMINAL 

************************************************************************ 


***  EXECUTION  CONTROL  STATEMENTS  *** 

********************  ************************************************ 

method  RECT 

timer  DELT=.?0,FI NT l Ms A9. 3, PR DEL =0.2, 3UTDEL=0.* 


*♦*  OUTPUT  CONTROL  STATEMENTS 


*** 


»RINT  OA(2-ll)tOBI2-ll) 

OUTPUT  I IME.CAI  UI.CBI  III 
PAGE  SR0UP=O.,?.E- 3) 

OUTPUT  TIME, OAT, OBT 
PAGE  GROUP=(D.,2.E-BI 
OUTPUT  T I ME , C A ( 1 1 ) , CB 1 1 1 ) 

PAGE  XV»L3T,HEIGHT*5.0,MlDTM«f».0,r,ROUP*IO.  ,2.E-3) 
OUTPUT  TIME, OAT, QBT 

PAGE  XYPL0T,HEIGHT=5.3,MlDTU=8.0.GR3UP=(0.,2.E-3) 


PNO 

STOP 
E NO JOB 


APPENDIX  D 


SINGLE -SOLUTE  P RE SATURATION  STUDIES 
SOLUTION  AND  SOLID-PHASE  STUDIES 


FIGURE  D-l : FLUID-PHASE  CONCENTRATION  PROFILE  FOR  DNOSBP 

FLU IDIZED-BED  ADSORBER. 


10.0 


FIGURE  D-2:  SOLID-PHASE  CONCENTRATION  PROFILE  FOR  DNOSBP  (79.56yM)  IN  A 

FLU IDIZED-BED  ADSORBER. 


IGURE  D-3:  FLUID-PHASE  CONCENTRATION  PROFILE  FOR  DNOSBP  (80.0yM)  IN  A 


(iuej6/sa  Loin  . 0lx)  uoi^ejq.u9ouoo  asend-P!-  L°S 


FUJI DIZED-BED  ADSORBER. 


■PHASE  CONCENTRATION  PROFILE  FOR  2,4-DCP  (80.0uM 
[ZED-BED  ADSORBER. 


mm 


r IGURE  D-7:  FLUID-PHASE  CONCENTRATION  PROFILE 

FLUIDIZED-BED  ADSORBER. 


FIGURE  D-8:  SOLID-PHASE  CONCENTRATION  PROFILE 

FLUIDIZED-BED  ADSORBER. 


FIGURE  D-13 : FLUID-PHASE  CONCENTRATION  PROFILE  FOR  2,4-DCP  (80.0uM)  IN  A 

PACKED-BED  ADSORBER. 


FIGURE  D-l 5 : FLUID-PHASE  CONCENTRATION  PROFILE  FOR  2,4-DCP  (80.0pM)  IN  A 

PACKED-BED  ADSORBER. 


.OyM) 


FIGURE  D-17:  FLUID-PHASE  CONCNETRATION  PROFILE  FOR  DNOSBP  (82. 

IN  A CSTR  ADSORBER. 


FIGURE  D - 1 9 : FLUID-PHASE  CONCENTRATION  PROFILE  FOR  DNOSBP  (80. 

IN  A CSTR  ADSORBER. 


FIGURE  D-21:  FLUID-PHASE  CONCENTRATION  FOR  2,4-CDP  (82.21yM) 

IN  A CSTR  ADSORBER. 


FIGURE  D-22:  SOLID-PHASE  CONCENTRATION  PROFILE  FOR  2,4-DCP  (82.21uM) 
IN  A CSTR  ADSORBER. 


FIGURE  D-23:  FLUID-PHASE  CONCENTRATION  PROFILE  FOR  2,4-DCP  (76.64uM) 

IN  A CSTR  ADSORBER. 


FIGURE  D-24 : SOLID-PHASE  CONCENTRATION  PROFILE  FOR  2,4-DCP  (76.64UM) 

IN  A CSTR  ADSORBER. 
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